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Abstract 

We consider particles on a one-dimensional lattice whose evolution is governed by 
nearest-neighbor interactions where particles that have reached size zero are removed 
from the system. Concentrating on configurations with infinitely many particles, we 
prove existence of solutions under a reasonable density assumption on the initial data 
and show that the vanishing of particles and the localized interactions can lead to non¬ 
uniqueness. Moreover, we provide a rigorous upper coarsening estimate and discuss 
generic statistical properties as well as some non-generic behavior of the evolution by 
means of heuristic arguments and numerical observations. 
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1 Introduction 

Discrete systems with local interactions are used in a wide range of models for growth 
processes. Examples are the coarsening of sand ripples [HK02] or the clustering in granular 
gases [vdMvdWLOl], where matter is transported between adjacent ripples or clusters, and 
certain particle hopping models [Spi70], where particles occupy cells on a one-dimensional 
lattice and jump between neighboring cells with rates depending on the particle density. 
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Figure 1.1: Numerical simulation of (1.1) with /3 = 0.5 and initial particle sizes uniformly dis¬ 
tributed in the interval (0,1). Left. Mean particle size (black) and maximum particle size (gray) 
over time, both in logarithmic scale. The thin line represents the power law t i—>■ Right. 

Particle size density scaled by the mean size at three different times during the evolution. 


Discrete systems with local interactions also appear as reduced models of more complex 
situations such as the coarsening of droplets in dewetting thin hlms [GW03, GW05], the evo¬ 
lution of phase domains in the convective Gahn-Hilliard equation [WORD03], or the motion 
of grain boundaries in polycrystalline materials [HNO04]. Furthermore, local interactions 
emerge in discrete or discretized ill-posed diffusion equations such as the Perona-Malik 
equation used in image segmentation [PM90] or population dynamics [HPO04], and there 
they lead to localized aggregation and coarsening as well [EG09, ES08]. 

In this paper, we consider a discrete system of particles of sizes (xj), j G Z whose 
evolution for time t > 0 is governed by 

Xj = - 2x7^ + xjf 1 (1.1) 

as long as xj >0; particles that have reached x = 0 are removed from the system and the 
remaining ones are relabeled. The parameter /3 is a positive real number. Equation (1.1) 
is closely related to the first examples above where our particles represent sand ripples, 
gas clusters or lattice cells. It can also be regarded as the discretized backward-parabolic 
equation Xj = Ag{xj) with decreasing mass transfer function g{s) = s~^ and A denoting 
the discrete Laplacian on Z. It differs from the discrete Perona-Malik equation, though, 
since the latter has a forward-parabolic region near the origin, whereby small particles do 
not vanish but cluster near the boundary to the backward-parabolic region. Moreover, 
due to its rather simple structure we see (1.1) as a useful toy problem for studying local 
correlations in coarsening systems and thus as a step towards understanding more complex 
and higher dimensional evolutions such as grain growth or the Mullins-Sekerka flow with 
positive volume fraction. 

The most interesting aspect of (1.1) and the examples above are the universal or generic 
statistical properties, which the evolutions seem to exhibit. Often, these can be derived 
formally from the equation or observed in numerical simulations of large finite systems, and 
for solutions to (1.1) the key findings are the following. 

1. Particles disappear in finite time, whereas the total particle size formally 

conserved during the evolution. Thus, the typical particle size grows, and a dimen¬ 
sional analysis suggests the growth law 

1 

typical particle size ~ tP+^. (1-2) 

The latter is confirmed by numerical results as depicted in Figure 1.1. 
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Figure 1.2: Top left. Sketch of initial data where one large particle, here xq = 10, is surrounded 
by small particles Xj = 1. Bottom row. Snapshots of the evolutions for /3 = 0.1 {left, black) and 
j3 = 0.3 {right, gray) at a later time. The small black particles vanish sequentially beginning at the 
large one, whereas the small gray particles become disordered. Top right. Maximum particle size 
for both systems over time, in logarithmic scale. For /3 = 0.1 the largest particle grows linearly in 
time (represented by the thin dotted line) instead of adhering to the generic growth law (1.2) (thin 
dashed line), while for fi = 0.3 its growth approaches the generic law (thin solid line). See Section 
2.3 for a more detailed discussion of the underlying dynamics. 

2. The simulations in Figure 1.1 also indicate self-similarity of the particle size distribu¬ 
tion in the sense 

particle size density(t, x) ~ t i^+^p(^xt 
for some profile p: [0,oo) —>■ [0, oo). 

3. Not all solutions display the behavior just described. For instance, constant initial 
data remain constant for all times, and as indicated in Figure 1.2 there are solutions 
where certain particles grow faster than the law (1.2). 

The analysis of coarsening in equations as (1.1) therefore comprises two problems: first, 
to derive an appropriate upper bound for the (typical) particle size that is true for all 
solutions, and second, to find conditions on the initial data that guarantee the generic 
behavior of corresponding solutions. The second issue is in general very hard to deal with 
and any mathematical consideration must be tailored to the specific evolution, whereas 
upper coarsening estimates have been obtained for several particle growth processes using 
numerical simulations, heuristic arguments or variants of the method introduced in [KO02]. 
Usually, however, it is difficult to apply these approaches rigorously to infinite particle 
systems. 

In the dynamics of (1.1), on the other hand, finitely many particles generically disappear 
in finite time, and due to vanishings and local interactions there is also no evident equation 
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Figure 1.3: Sketch of the initial data for our non-uniqueness result. Large particles x^j ^ 1, j > 0 
are followed by small ones R-ij+i ~ R‘ij +2 Rj where —>■ 0 as j —>■ oo. The small particles 

vanish pairwise from right to left, and thereby tiny perturbations of a large particle (thin gray error 
bars) are transported to the left and amplified. Non-uniqueness is then a consequence of the infinite 
number of particles, which allows arbitrarily large distances to be covered and thus amplifications 
of order 1. 


for the evolution of a particle size distribution. In order to rigorously characterize the 
statistical properties, our goal is thus the mathematical analysis of infinite systems. In 
the current paper, we study well-posedness of (1.1) for infinitely many particles and prove 
existence of solutions under a density condition on the initial distribution. More precisely, 
assuming that there are two constants L > 0 and d > 0 such that 

> d 

jeR 

for any region R <Z "L containing L consecutive particles we show that there are locally 
Holder continuous trajectories Xj-. [0,oo) —)■ [0,oo) that satisfy (1.1) in an integral sense 
and almost everywhere as long as Xj > 0. In particular, a solution contains infinitely many 
positive particles at any time t > 0. 

It has been observed that infinite particle systems might have pathological solutions, for 
instance in [Lan75] where the author discusses a Hamiltonian system of classical particles 
interacting by means of hard core potentials. He demonstrates that for well-prepared initial 
data it is possible to transfer energy from infinity to a finite region in space via infinitely 
many particle collisions in finite time, and along a similar idea we show that solutions 
to (1.1) are not uniquely determined by their initial data. To this end, we consider a 
conhguration as sketched in Figure 1.3 where large particles x^j ~ 1 are followed by two 
small ones x^j+i ~ X 3 j +2 ~ Rj with —>■ 0 as j —>■ oo. Exploiting the effects of vanishings 
and the localized interactions we prove that (even in finite systems) tiny perturbations 
pass through the conhguration and are amplihed during the evolution. Arbitrarily small 
initial perturbations may reach a size of order 1 if the distance they travel, and thus their 
total amplihcation, is sufficiently large. In the inhnite system, we then obtain two different 
solutions by considering unperturbed data as well as decreasing perturbations at increasing 
particle indices. 

As a hrst step towards understanding the statistical properties of (1.1), we provide a 
simple upper coarsening estimate relating the mean particle size of a system and the time 
that is necessary to transport a fraction of it from vanishing particles to remaining ones. 
We also discuss part of the non-generic behavior indicated in Figure 1.2 relying on heuristic 
arguments for the mass transfer in the system and numerical investigations. 
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The rest of the paper is organized as follows. In Section 2 we introduce our precise 
setting, which in particular comprises how we deal with vanishing particles. We also collect 
our coarsening results in Sections 2.2 and 2.3. The existence theorem and its proof are 
contained in Section 3, while Section 4 is devoted to non-uniqueness. 

2 Setup and coarsening 

A precise description of the dynamics of (1.1) requires that we keep track of particles 
and their neighborhoods throughout the evolution. For this reason we consider sequences 
X = (xj) in the state space 

= {x = (xj) : 0 < Xj < C for all j G Z and some constant C > 0} , 

which include vanished particles Xj = 0 as well as living ones Xj >0; the term particle can 
refer to the size Xj as well as the index j. Since only adjacent living particles interact with 
each other, we use the term neighborhood of a particle Xj to denote the two nearest living 
particles where 

(T_(j,x) = sup{/c < j : Xfc > 0} , a+(j,x) = inf {/o j : Xfc > 0} ; (2.1) 

in the case that one of the sets in (2.1) is empty we set sup 0 = —oo, inf 0 = oo and x±oo = 0. 
We also write 


cj_(oo, x) = sup {A: : Xfc > 0} , (T+(—oo, x) = inf {/c : x^ > 0} 

for the largest and smallest particle index in a finite system, respectively. 

Given a function g: [0,oo) —)• [0,oo), we denote by 

^crg{Xj) j.j^(Xj) (ji a - (j ,x)) 2gf(Xj) -|- (7(^(T+(j,a:))^ ' X{xj>0} 

the living-particles-Laplacian of the sequence g{x)\ for simplicity of notation we omit the 
indices j and x in cr, a± when the meaning is clear from the context. In our model problem 
we have 

g{s) = s~^ for 0 < s < oo, ^(0) = 0, 
where /? > 0, and with the convention 0“^ = 0 we write (1.1) as 

Xj{t) = A(^g{xj{t)) = AfjXj{t)~^ for j G Z, t > 0 (2.2) 

or 

x{t) = Aa-x(t)“^ for t > 0. 

By definition of A^, the right hand side of (2.2) is well-defined for all j G Z. In particular, 
we have Xj{t) = 0 if Xj{t) = 0 and Xj{t) = Xj-i{t)~^ — 2xj{t)~^ if Xj-i{t) > 0, Xj{t) > 0, 
but a+{j,x{t)) = oo. Finally, if Xj{t) > 0 for all j G Z, then A^ is the standard discrete 
Laplace operator on Z. 

2.1 Notion of solution 

We now define our notion of solution. As we will prove existence by passing to the limit in 
the number of particles in a system, an integral formulation is appropriate. On the other 
hand, for a priori estimates the differential equation (2.2) is more convenient and hence we 
briefly discuss equivalence of both formulations. 
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Definition 2.1 (Integral solution). For T G (0, oo] we call x = {xj)j^z- [0,F) —a 
solution to (2.2) with initial data x(0) if 

1. each Xj is continuous in [0,T); 

2. each function xJ^X{xj>o}: where 0“^ • 0 = 0, is locally integrable in [0,T); 

3. for any j G Z the equation 

rh 

Xj{t2) - Xj{ti) = / A„Xj{s)~^ds (2.3) 

Jti 

is satisfied for all 0 < ti < t 2 < T. 

For a solution x as in Definition 2.1 and m, n G Z, m < n, we denote by 

n 

Mm,nit) = ^ Xkit) (2.4) 

k=m 

the mass of the particles m,... ,n at time t G [0,T). Using (2.3) and the definition of Aq-, 
we compute 

Mm,n{t2) - Mm,nitl) = f “ ^a+{a.{ni)) ~ ^C7^(c7+{n)) + ^(7+{n)} ' X{M>0} ds (2.5) 

for 0 < ti < t 2 < T, where Mm,nit) > 0 implies that there is at least one living particle 
among m,... ,n and thus that fT_(m) <m < iT+(cr_(m)) < (T_(cJ+(n)) <n< cr+(n) in the 
integrand on the right hand side. 

The models in the introduction assume that vanished particles do not reappear at later 
times, and with the help of (2.5) we can show that our solutions indeed have this property. 

Lemma 2.2. Ifxjiti) = 0 for some j G Z and ti G [0,r), then xjit) = 0 for all t G [ti,T). 

Proof. Assume for contradiction that the claim is wrong and denote by 

t 2 = inf {t > ti : Xjit) > 0} G [ti,T) 

the time up to which Xj remains vanished. Suppose first that m = cj_(j, x(t 2 )) and n = 
cr+(j, xit 2 )) are both finite and fix ta G it 2 ,T) so that by continuity of particle sizes we have 

Xmit)~^ < c, Xnit)~^ < c and Xkit)~^ > 2c, fe = m + 1,..., n — 1 

for all t G [^ 2 ,^ 3 ] and some constant c > 0. The change of Mm+i,n-iit) then satisfies 

Xmit)~^ - X„^^m,x{t))it)~^ - X„_^ri,x{t))it)~^ + Xnit)~^ < -2c, ( 2 . 6 ) 

for t G [^ 2 ,^ 3 ] such that Mm+i,n-iit) > 0, hence (2.5) gives 

Mm+l,n-lit) < -2c|{s G it 2 ,t) : Mm+l,n-lis) > 0}| 

for all t G [t 2 ,t 3 ]. This is a contradiction to 0 < Xjit) < Mm+i,n-iit) for t sufficiently close 
to t 2 . 

In the cases that m = —oo or n = oo or both, the corresponding positive term(s) in 
(2.6) disappear(s). Since, however, the negative contribution —x„^(^m)it)~^ — x^_(^m)it)~^ 
remains, we still obtain the contradiction as above. □ 
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Remark (Vanishing times). As a consequence of Lemma 2.2 each particle Xj has a unique 
vanishing time tj G [0, T) U {oo} such that 

Xj (t) > 0 for t < Tj and Xj (t) = 0 for t > Tj , 

where we set tj = oo if and only if Xj{t) > 0 for all t G [0,T). Moreover, the functions 
t I—7- cr±{j,x{t)) are monotone, right-continuous and change their value in vanishing times 
only. 

Lemma 2.3 (Solving the ODE). For j G Z let Vj = {t : a±{j,x{t)) is not continuous} be 
the set of times when a neighbor of Xj vanishes. Then Xj is continuously differentiable for 
all t G (0, T) \ ({tj} U Vj) and satisfies (2.2). 

Proof. Fix j G Z. For t > Tj the claim is trivial and for t G (0,rj) \ Vj it follows from 
the Fundamental Theorem of Calculus applied to (2.3) in a sufficiently small time interval 
around t. More precisely, with m = a-{j,x{t)) and n = a+{j,x{f)) hxed, continuity of Xj, 
Xm and Xn provides <5 > 0 and c > 0 such that 


min{xm('S), a:j(s), Xn(s) '. s G [t — 5,t + (5)} > c > 0. 


Hence, (T_(j, x(s)) > m and cj+(j, x(s)) < n for all s G {t — 6,t + 6), which implies that as 
functions of s both change their value only at finitely many times in Vj D (t — 5,t + 6). As 
t is not such a time, we may decrease 6 to obtain a-{j,x{s)) = m and (T+(j, x(s)) = n for 
all s G (f — (5, t -b 5). Therefore, 


Xj{t + h) — Xjif) 

h 



2Xj{s) ^ + Xn{s) ^ds 


for all \h\ < (5, and the claim follows as stated above. 


□ 


Remark. The preceding lemmas show that every solution x has the property that 


1. for each j there exists a unique vanishing time Tj ; 

2. each Xj is continuously differentiable in (0,r) \ V, where V = {xk : k G Z}, and equa¬ 
tion (2.2) holds. 

The other way round, if x : [0, T) —>■ Iff satishes the second item above with V replaced by 
an arbitrary null set of [0,T) and Definition 2.1(1 and 2), then also the integral identity 
(2.3) is true. Moreover, uniqueness of vanishing times follows directly from the differential 
equation if the latter holds almost everywhere, since the argument in the proof of Lemma 
2.2 is easily adapted. 


2.2 Upper coarsening estimate 

Equation (2.2) can formally be interpreted as the H~^ gradient flow on living particles of 
the energy 

E{x) = X] if ^ / 1 or E{x) = ^ Inx^ if ^ = 1, 

k\ Xk>0 k: Xk>0 

and such a structure has proved successful for obtaining upper coarsening bounds in similar 
(finite) systems; see for instance [EG09, ES08] for the application of [KO02] to a discrete 
H~^ gradient flow. However, besides the fact that energy and metric tensor of infinite 
systems are in general inhnite, this gradient flow structure does not seem useful here, since 
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(even in finite systems) at each vanishing time the local correlations in the metric tensor 
change and for (5 > 1 the energy tends to negative infinity. We thus take a different 
approach, which regards an upper coarsening estimate as a lower bound on the time needed 
to transport a certain amount of mass through the system and which is derived from the 
basic inequality Xj > —2x~^. 

For illustration, suppose that x = {xk)k=-n,...,n is a solution to (2.2) with 2n + l nonzero 
initial particles whose mass is M = X]|fc|<n®fc(0) denote by 


T = inf < t > 0 : ^^^(O) > -M 

[ Xk{t)=0 

the smallest time such that half of the initial mass has been transported from vanished 
particles to living ones. Then, there must exist a particle Xj with vanishing time Tj < T and 
initial size at least M/(2(2n + l)), because otherwise we would have Ylx^(T)=o^kiO) < M/2. 
From Xj > —2xJ^ we infer that Xj{0)^^^ < 2{f3 + l)rj and hence obtain the estimate 

< x,(0) < (2(/3 + 1 )t,)^ < [2(/3 + l)]im r?Ti, 

which bounds below by the initial mean particle size and a constant depending 

on /3 and the fraction of mass that is transported. 

Assuming that the mean particle size initially exists, we prove a similar estimate for 
infinite particle systems. 


Proposition 2.4. Letx: [0,cx)) —)• be a solution to (2.2) which initially has finite mean 

particle size 

1 

X = lirn Xn G (0,oo), Xn = ^ ^ 3:^(0) 


and let 


as well as 


2 n + 1 


j=-n 


T„, = inf 


t > 0 : 


1 


I ^i(O) — 2^"' 


2 n + 1 

j=—n,...,n 
Xj(t)=0 


r = inf < 


t > 0 : liminf 


1 


n^oo 2n + 1 


Xj{0) > ^x 


j=—n,...,n 

Xj(t)=0 


. 


Then there is a constant c > 0 such that lim inf,i_>.oo Tn > and T > 


In Proposition 2.4 both, T and liminf are a measure for the time needed to transport 
half of the initial mean particle size from vanishing particles to living ones, but we do not 
know which one is better or more useful. Furthermore, the result only applies when suffi¬ 
ciently many particles do indeed vanish. We expect that this is the generic case for long-time 
evolutions, since the discrete backward-parabolic equation (2.2) increases differences among 
particle sizes and since small particles surrounded by larger ones disappear sufficiently fast; 
see Lemma 3.6. On the other hand, there are situations such that no or very few particles 
vanish and the possible rate of coarsening is not captured by our estimate; consider for 
instance a system of large regions of equal particles and compare the mass concentration 
example in Section 2.3. 
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Proof of Proposition 2.4- As for finite systems above we obtain Tn > cXn'^^, and taking 
the lower limit as n —>■ oo proves the first inequality. 

By definition of T there is for any e > 0 a time G [T, T + s] such that 


hminfK(i£) > 7 ;^, 

n^oo Z 


ji=—n,...,n 
Xj{t)=0 


Hence, there is no = no(e) such that Vnite) > X/2 — e for all n > no and the argument for 
hnite systems yields 

1 1 

-X -e< < c{T + e) /s+i. 

Sending e —)• 0 finishes the proof. □ 


2.3 Example: mass concentration versus mass spreading 

We conclude this section with a heuristic discussion of the example shown in Figure 1.2 
where initially a large xo(0) 3> 1 is surrounded by Xj(0) = 1, j / 0. 

The particle xq clearly grows during the evolution, and assuming that it is sufficiently 
large we may neglect its contribution to the particle velocities. Furthermore, due to the 
symmetry x-j = xj for j G Z it suffices to consider only the particles xj, j > 0, and the 
dynamics of the example are effectively governed by 

xo ~ xi ~ -2xf^ + xf^, Xj ~ xj^^ - 2xJ^ + xj^^ for j > 1. (2.7) 

We let 

Xj{t) = 1 + s = /3t 

and linearize 

xj^{t) = exp{—pinXj{t)) ~ 

so that (2.7) becomes 

Ano~2e-“S ^ _2e-“i + e-“^ ^ - 2e-“^'+ (2.8) 

with initial data uo{0) 1 and Uj{0) = 0, j > 0. The initial conditions immediately show 

that 

A«o(0)~2>0, A^^(o)^_i<o, £i^i(0)~0 for j > 1, 

and by induction it is easy to see that 

~ 1 > 0’ (s)^A(0)=0 k>j>l. 

Indeed, for j > 0 equation (2.8) reads ^ ~ where A is the discrete Laplacian on Z, 

and taking the A:-th time derivative we obtain 






which enables the induction argument. As a consequence, we find (—l)-^Uj(s) > 0 for small 
positive times s, and using these inequalities in (2.8), we conclude that they remain true 
for all s. Moreover, a comparison argument for Uj +2 — Uj shows that uq > U 2 > u^ - ■ ■ and 
■ ■ ■ P U 5 P U 3 P ui- In terms of the particles Xj this mean that those with an even index j 
grow while those with an odd index shrink as long as (2.8) is a valid approximation of the 
original dynamics. 
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We now consider two cases, namely /? being very small and very large, respectively. For 
small /3, the requirement Uj = 0(/3) for j > 1, which means that none of the corresponding 
particles has vanished, and (2.8) yield 

£uo ~ 2 and £ui ~ -1. 

The key findings are the following. 

1. The particles xi and x_i vanish when ui reaches a value of order —/?, and that happens 
at a time of order s ~ /3 or f ~ 1. 

2. On this time scale the particles Xj, |j| > 1 hardly change size due to \uj\ = 0(/3). 
Therefore, after xi and x_i have vanished, we are in the same situation as initially 
where xq is surrounded by small particles Xj « 1, now indexed by \j\ > 1. 

3. The mass xi(0) +x_i(0) = 2 is almost completely transported to the particle xq, and 
the latter grows linearly in the time t. 

With regard to the second and third item we refer to the case of small f3 as sequential 
vanishing or mass concentration. 

In the second case, when /? is large, we use that —ui has to become large in order that 
xi vanishes. Due to U 2 > 0, we have 

~-2e-“i + <-2e-“S 

which implies that the time for ui to reach a size of order —/3 is of order s ~ 1. Moreover, 
we have 




- 2 e "“2 + e 




-U3 > _ 3 

^ ~ 2 ds “1 9 


due to ^2 > 0 > rt 3 , and we conclude that 

Ms) > -iMs) - is- 


(2.9) 


The main observations are now the following. 

1. The particles xi and x_i vanish at a time of order s ~ 1 or t ~ 

2. The change of the other particles is not known to be smaller. In fact, by (2.9) the 
change of U 2 is of the same order as the change of ui ~ —/? at times of order s ~ 1. 

3. Due to (2.9) and the equation for uq in (2.8), the mass xi(0) = 1 is split between xq 
and Xj, j > 1 with at least half of it (up to an error 1/(3) going to Xj, j > 1. Similarly, 
the mass of x_i is split between xq and Xj, j < —1. 

In view of the third item we call this case mass spreading. 

Our heuristic observations are confirmed by numerical simulations. In fact, the latter 
even yield a refined picture which suggests that there is a critical value (3c ~ 0.27 such that 
for initial data xq = Xq > 1 and Xj = 1, |j| > 1 one of the following three cases holds; 
compare Figure 2.1. 


1. If (3 > Pc then for all Xq > 1 mass spreading occurs. 

2. For each P < Pc there is an Xg such that for Xq < Xg mass spreading takes place. 

3. For P < Pc and Xq > Xp mass concentration arises. 

Moreover, our numerical results also indicate that in the mass spreading case the particle 
sizes become disordered over time as indicated in Figure 1.2. Currently, however, we are 
not able to provide good heuristics for or a suitable notion of being disordered. 
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p 

Ao 


mass 

mass 


spreading 

concentration 

.050 

1.72 

1.73 

.100 

8.00 

8.06 

.125 

28.00 

28.20 

.150 

84.60 

85.40 

.175 

188.00 

189.00 

.200 

537.00 

541.00 

.225 

2456.00 

2470.00 



Figure 2.1: “Phase diagram” for the example in Section 2.3. There is a critical /3c, which according 
to simulations lies between 0.26 and 0.27, such that for /3 > /3c always mass spreading occurs. If 
/3 < /3c, on the other hand, we have mass spreading for small Xq and mass concentration for large 
Xq. Some numerical values are shown in the table and as gray dots in the diagram. 


3 Existence of solutions for initial data with positive density 


In this section we prove existence of solutions for suitably distributed initial data which 
contain no large regions in Z of vanished or very small particles. We also suppose that 
initially there are no vanished particles. 

Assumption 3.1. All particles in the initial conhguration x = (xj) G are positive. 
Moreover, there exist two constants d > 0 and L > 0 such that Z can be subdivided into 
regions of length L with mass density at least d; that is, there is /cq G {0,1,..., L — 1} such 
that for 

R{k) := {ko + kL,... ,ko + {k + 1)L — 1} , /c G Z 

we have 

^ Xj> d. (3.1) 

jeR{k) 

Without loss of generality we assume ko = 0. 

Remark (Traps). Equation (3.1) implies the existence of indices jk G R{k) such that Xj,, > d 
for all G Z. In particular, we have jk —t ±oo as —>• ±oo and 0 < jk+i — jk < 2L. We 
call the Xj^ and jk traps, and given j G Z we denote by j± the two nearest traps (unequal 
to j itself), that is j_ = max {jk < j} and j+ = min {jk > j}', obviously, we have 

j+ - j- < 4L, j+ - j < 2L, and j - j_ < 2L 


for all j G Z. Finally, if x is a particle conhguration with traps (jk) as above, then there 
is at least one trap among 2L consecutive particles. Thus, the existence of traps is up to a 
factor for L equivalent to the density inequality (3.1). 

An example for a particle ensemble that satishes Assumption 3.1 is a positive x G 
such that the one-sided mean particle sizes 


lim - 

n—>-oo n -|- 1 


n 


j=0 


lim -- Xj 

n^oo n + 1 

j=-n 
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exist and are positive. To see this, assume for contradiction that there are no appropriate 
traps, that is, for any 6 > 0 the indices (jk) of particles Xj,, > 5 either are bounded from 
above or below or have unbounded distance jk+i — jk as k ^ +oo or /c —>■ —oo. If the 
latter happens as k ^ oo, we can hnd an index /cq = ^o(<^) such that jk+i — jk > 1/^ 
for all k > ko, which implies # {jk ■ jko ^ jk <n}<{n + l)(i for n > jk^- Together with 
# {jk : 0 < jfc < jko) < {jko + 1) and xj < C for all j G Z we obtain 


n + 1 


j=0 


n + 1 


E 


Xj + 


.j=0,...,n j=0,...,jko 


Xj -\- 


i>S 


E ^ 

Xj>5 


<6 + ^ + C8 

n + 1 


for n > jko, and sending n —)• oo gives a contradiction as (5 > 0 is arbitrary. The other cases 
are similar. 

The sequential vanishing example in Section 2.3 provides strong evidence that, at least 
for small {3, a condition like Assumption 3.1 is necessary for bounded solutions to exist. 
Consider an initial configuration that contains infinitely many regions Rj where a large 
particle of size 1 is surrounded by Lj particles of size Sj on each side, and assuming that x is 
a solution for such initial data, rescale time and particle size by means of y{t) = x{£j~^^t)/ej. 
Then, y still satisfies equation (2.2), particles of size x = £j are transformed to y = 1 and 
those of size x = 1 to y = 1/£j. The particles y in Rj thus resemble exactly the sequential 
vanishing situation in Section 2.3, provided that /3 and £j are small and Lj is sufficiently 
large so that the particles around the large one in Rj are not influenced by the adjacent 
regions Rj±i. In particular, the large particle in Rj gains all the mass of the nearby vanishing 
ones and its size y grows linearly in time. In terms of x, the large particle grows like 


.-(/ 3 + 1 ) 


t£j 


= tej^, 


which means that we can arrange an increase of order Ej in a time of order if 

LjEj S> Ej . Hence, by Ej —)• 0 and choosing Lj appropriately we obtain an initial 
conhguration x G so that some particles become arbitrarily large in arbitrarily small 
time. 


3.1 Local-in-time a priori estimates 

We now derive key properties of solutions with initial data as in Assumption 3.1. In the 
following, we let x: [0,T) —)■ be such a solution as in Definition 2.1 and for simplicity 
of notation we consider times 0 < t < T without indicating it explicitly. Throughout the 
paper, we denote by C and c generic constants that depend on /3 only. 

According to Lemma 2.3 the differential equation (2.2) is true for all but the vanishing 
times V = {r^ : A: G Z}, and dropping the non-negative contributions of Aq- gives 

Xj{t) > —2xj{t)~^ 

whenever Xj > 0 and f 0 V. Integrating from ti to t 2 , where 0 < < t 2 < Tj, we find 

Xj{t2f+^ - Xj{tif+^ > -2{/3 + l){t2 - h) (3.2) 


and conclude that 


Xj { t 2 ) > (xj(ti)^+^ - 2 { f 3 + l)(t2 - ti 


1 

^ + 1 


(3.3) 
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as long as the difference on the right hand side is non-negative. Hence, we obtain the 
half-life estimate 

Xj{t) > for all t <ti Cxj{ti)f^~^^ 

and in particular the persistence estimate for traps 

Xjkit) > - for all t <T* := ,T). 

The latter implies 

(^-U,xit)) > j-, a+{j,x{t)) < j+ for all t < T*, j € Z, 

and therefore the number of different particles that up to time T* can appear as neighbors 
of any hxed particle xj is at most 4L. 

Next, we consider the local masses Mm,nit) = J2k=m^kit), m < n as dehned in (2.4). 
Similar to (2.5), we have 

^^Mm,nit) = X„_l^m){t)~^ “ X„_(„^(m))it)~^ “ 2:^+(<7_(n))(i)"^ + X„^(ri)it)~^ (3.4) 

whenever Mm,nit) > 0 and f 0 V, and we infer that 

^^Mm,nit) < X„_(m)it)~^ + X^^(^n)it)~^ “ 2Mm,nit)~^ (3.5) 

since 0 < x^_(^_^(m))(t) < Mm,nit) and 0 < X(^_^(^_(„))(f) < Mm,nit). Furthermore, if jk < ji 
are two traps, the persistence estimate and integration of (3.4) provide 

fd\~^ 

^jkJii^'^) — ^jkJii^^) ~ ^ (^ 2 / ~ (^•®) 

and 

fd\~^ 

^ife+iji-i(^2) < Mj,^+ij,_i(fi) -h 2 y-j it2 - h) (3.7) 

for all 0 < ti < t 2 < T*. From the inequalities (3.6) and (3.7) we draw several conclusions 
that are essential ingredients of our existence result. First, the lower bound (3.6) implies 
that the mass density inequality (3.1) in Assumption 3.1 is stable up to time T* in the 
following sense. 

Lemma 3.2. For any q € N, q > 1 and any k € Z the mass density in the region 


Rqik) := {qkL, ..., qik -|- 1)L — 1} = Riqk) U • • • U Riqik -|- 1) — 1) 


of length qL satisfies 


for all t e [0, T*]. 


iz ^ 

j£R^{k) ^ 


1 + C/L \ 
Q ) 


d 


Proof. Each of the q regions Riqk), ..., Riqik -|- 1) — 1) of Assumption 3.1 contains a trap, 
and denoting by G Riqik -|- 1) — 1) and k~ G Riqk) the largest and smallest of these we 
infer from (3.6) that 


^ Xjit) > Mk-,k+i0) -Cd ^t> Mk-,k+i^) - Cd 

j£Rq[k) 
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for all t < T* < Moreover, Assumption 3.1 also implies 

q{k+l )-2 

E E Xj > (q — l)Ld, 

n=qk+l jeR{n) 

and combining both inequalities finishes the proof. □ 

The key point of Lemma 3.2 is the trade-off between the lower bound and the number 
of regions in the density estimate: by choosing q sufficiently large, we can ensure that the 
initial bound d is almost preserved. 

From the upper mass estimate (3.7) we derive Holder continuity and boundedness of 
the particle trajectories. 

Lemma 3.3. We have 

\xjit 2 ) - Xj{ti)\ < C{L + l)\t 2 - ti\^ 

for any j G Z and all Q <ti <t 2 <T*. 

Proof. Given j G Z, it suffices to consider t 2 < T* such that Xj{t 2 ) > 0. Using (3.7) with 
the two nearest traps j± and T* < Cd^~^^ we have 

A/j-_+ij+-i(t2) < + 2{d/2)~^T* P+^{t2 — ti)'3+i 

< Mj_+ij+_i(ti) + G(t2 -H)^, 

and separating Xj from both mass terms yields 

Xj{t 2 ) - Xj{tl) < ^ [Xk{tl)-Xk{t 2 )) +C[t 2 -ti)l^. (3.8) 

j-<k<j+ 

We split the sum in (3.8) into large and small particles, respectively, and consider the 
summation over 


UL = {k£ {j- + {j} : Xkfti) > s} , 

Us = {kG {j- + l,...,j+-l}\ {j} : Xk(ti) < e} , 

separately, where e is defined by = 4(/? + l)(t 2 ~ ^i)- For the small particles we simply 
estimate 

Xk{ti) - Xk{t 2 ) < e{j+ - j-) < CL{t 2 - (3.9) 

keUs 

using Xk{t 2 ) > 0, — j- < 2L and the definition of s. For k G Ul, on the other hand, we 

hnd Xk{ti)^~^^ — 2(/3-|-l)(t 2 —^i) > 0, hence we may use (3.3) and the elementary inequality 
— 6^+^) > (a — 6)^+^ for /? > 0 and 0 < 6 < a to get 

Xk{t2) > Xkih) - C{t2 - ti)U^ 


and conclude 

Xkih) - Xk{t 2 ) < CL{t 2 - (3.10) 

keUL 

Combination of (3.8)-(3.10) yields the Holder estimate in the case Xj(t2) > Xj{ti), while 
for Xj{t 2 ) < Xj{ti) the result is an immediate consequence of (3.2) and once more the 
elementary inequality > (a — 6)^+^. □ 
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An upper bound on the particle sizes now follows from either the mass estimate (3.7) 
or the Holder inequality. 

Corollary 3.4. We have 


Xj{t) < + Cd and Xj{t) < Xj{0) + CLd 

for all t G [0, T*] and any j G Z. 

Finally, we characterize the behavior of vanishing particles. Inequality (3.2) is in case 
Tj < oo also true for t 2 > tj and immediately provides 

Xj{t) < C{Tj — t)/3+i for all 0 <t < Tj, (3-11) 

which can be understood in two ways. On the one hand, (3.11) is an upper bound on the 
vanishing particle Xj{t) in terms of Tj — t, and on the other hand, setting t = 0 gives a lower 
bound on the vanishing time Tj in terms of Xj(0). In the following two lemmas we prove the 
opposite assertion by showing first that the size of a vanishing particle is estimated from 
below by a power law and second that sufficiently small particles do indeed vanish. 

Lemma 3.5. There is a constant K = K{l3,d,L) such that 

Xj{t) > K{Tj-t)W^ 

for all t G [0, Tj] and any j G Z with Tj < T*. 

Proof. We prove the claim by contradiction. If it is wrong there exist solutions x” : [0, T"") —?■ 
n G N, such that each has a particle j”’ with vanishing time T^n < min(r’^, Cd^~^^) and 
a bad time < T^n with 


(r- 


as n —)• oo. 


We may shift the particle index of x” so that 0 = j” < j"' < < 4:L and then select 

a subsequence of n —)• oo in which and jlf are constant. In the following, we restrict 
our attention to this subsequence of n and to the particles fc = 0,..., j”; for simplicity of 
notation we drop the index n in the notation. 

Rescaling time and size we consider 


Vkis) 


Xkjt) 

{Tj - 4 ) ^ 


t = Tj + {Tj - tif)s, 


whereby the bad times 4 are translated into s = — 1 and the vanishing times of Xk into 
vanishing times of yk', in particular, we have aj = 0 and deduce 


yj{s) > 0 for all s G [—1, 0), yj{—l) —)• 0 as ra —>■ oo. 

It is clear that ^yk{s) = Ao-yA:(s)~^ whenever Xk{t) is differentiable at the time t corre¬ 
sponding to s. Moreover, since Xk is Holder continuous we also obtain 

\yk{s2) - yk{si)\ < C\s2 - 


for all si, S 2 G [—1,0] and k = 0 ,..., where C depends on (3 and L. 
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To achieve a contradiction we study the rescaled mass in a neighborhood of j in which 
the particles vanish as fast as Xj. For that purpose, we let B C {1,..., — 1} be the 

largest set of consecutive indices containing j such that 


liminf yk{ — l) = 0 for all k £ B] 

n—¥oo 


writing B = {m + 1,... ,l — 1} and restricting ourselves to a subsequence of n we may 
assume that the lower limits are attained as limits and that > c, y/(—1) > c for all 

n and some constant c > 0. Holder continuity then implies 

Vkis) > c — ^ for A: = m, I, 

s G [—1, —1 + (5) and sufficiently small <5 > 0, whereas the rescaled mass M(s) = Vkis) 

satisfies 

M(s) < M(-l) + 4LC6^ < 2"? I 

for s G [-1,-1 + 5), sufficiently small 5 G (0,1) and all sufficiently large n. Choosing 
6 = 6{f3, L, c) and n = n{f3, L, c, 5) so that these two inequalities are true, we obtain 

£m{s)= ^aVkis) <ymis)~^+ yi{s)~^- 2 M{s)~<^ < 

k&B 

yk(s)>o 

for all s G (—1, —1 + 5) \ {a^ '■ k G B} and therefore 

0 < < M(-l)^+^ - (/3 + 1)(5. 

The contradiction now follows from increasing n until M(—1)^^^ < (/3 + 1)5. □ 

Lemma 3.6. For any to G (0,T*) there is a constant rj = r]{(3, d, L,to) such that Xj{t) < y 
for some t G [0,to] <^'nd j G Z implies Tj <t + (xj(t)/iF)^"’“^, where K = K{f3,d,L) is the 
constant from Lemma 3.5. 

Proof. Lemma 3.5 yields tj <t + {xj{t)/K)l^~^^ provided that tj < T*; hence, we only have 
to prove the latter inequality in the case Xj{t) > 0. Furthermore, by a shift in time it is 
sufficient to consider t = 0. 

Following the same idea as in the proof of Lemma 3.5, we assume for contradiction that 
the claim is wrong. Then there are a sequence 0 as n —)> oo and solutions x^: [0, T") 

such that T” > to and 

a^j„(0) < rjn, > 0 for all t < min(T”, Cd^~^^) 


for some particle j”. Again, by shifting the particle index we select a subsequence of n —)• oo 
such that jT =0, j” and j” are constant, drop the index n in the notation, and restrict 
our attention to the particles A: = 0,..., j+. 

We now rescale time and particle sizes according to 


yk{s) 


Xkjt) 

V 


t = y^^^s 


and denote by H = {m + 1,...,/ — 1} C {0,..., j+} the largest set of consecutive particles 
containing j such that 


liminf?/fc(0) = Ok £ [0,oo) for all k £ B] 

n^oo 
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restriction to a subsequence of n allows us to assume that the lower limits are attained as 
limits and that for k = m, I we have yfc(O) —)> oo as n —>■ oo. As before, Holder continuity 
implies 

1 

yk{s) > c — for k = m,l, c = min ?/a;(0), 

k=m,l 

and s G [0,5) where 5 = while the rescaled mass M{s) = YlkeB Vkis) satisfies 

M{s) < M{0)+ALCs^. 

Since 5 —)■ oo, c —>• oo and M(0) —)• YlkeB oik SuS n ^ oo, we can for any 5 > 0 find an n 
such that 6 < S and 

Vkis) >c- C6^ > 2^/^(M(0) + 4LC6^) > 

for k = m,l and all s G [0, <5]. Hence, we obtain 

iM{s) < ym{s)-^ + yi{s)-^ - 2M{s)-^ < -M{s)-^ 

for all s G (0, 5) \ : /c G H} and 

0 < yj{ 6 f^^ < < M(0)^+^ - (/3 + 1)A 

Choosing 5 and n sufficiently large gives a contradiction. □ 

Remark. If Tj > T*, then Xj{t) > c > 0 for all t G [0,T*]. This and Lemma 3.5 if Tj < T* 
imply that every ^X{xj>o\ is integrable in [0, T*]. Hence, for initial data as in Assumption 

3.1 any solution in [0, T*) to the differential equation (2.2) as in the remark after Lemma 2.3 
is a solution according to Definition 2.1. 

3.2 Existence of solutions 

We now prove existence of a solution to (2.2) by truncating the initial data at a particle 
index n and sending n —>■ oo. To this end, given x°° G as in Assumption 3.1, we first 
look for a solution of the finite dimensional initial value problem 

x'j{t) = AcrX^{t)~^ if t > 0 and x^{0) = x^ for j = —n,...,n. (3.12) 

Recall that (3.12) is well-defined by definition of Aa-, in which contributions from beyond 
the two outermost living particles are dropped. 

By a solution to (3.12) we mean a continuous x”: [0, oo) —)■ A”, where 

= |x = (xj) G : 0 < Xj < C for j = —n, ... , n and some constant C > O} , 

such that x” attains the initial data in (3.12) and has the following property: For each 
j = —n, ... , n there exists a vanishing time G (0, oo] as defined before, and between 
t = 0, consecutive vanishing times and t = oo each x” is continuously differentiable and 
satisfies the differential equation x” = Ao-(x”)“^. In short, x^ solves the same differential 
equation as in the full problem, but with truncated initial data, and consequently, the 
results of Section 3.1 are true for x” provided that n is larger than the particles involved. 
In particular, we have 

1. persistence of traps 

^ 

for all t G [0,T*] and jk,n such that n > \jk\', 
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2. the upper bound 

x]{t) < C 

for t £ [0, T*] and all j, n such that n > \j± \ where C depends on I3,d, L and sup^g^ 1 

3. the Holder estimate 

\Xj(t2) - X'j{ti)\ < C\t2 - ti\^ 

for all 0 < ti < t 2 <T* and all j, n such that n > |j±|, where C = C{I3, L); and 

4. the vanishing estimate 

x]{t) < r] for some t < to xj <t + < T* 

for to £ (0,T*), sufficiently small rj = T]{j3,d, L,to) and all j,n such that n > |j±|, 
which in turn implies 

x]{t) > C{t^ - t)i^ 
for t G [0, r”], where C = C{f3, d, L). 

Furthermore, equation (3.5) for the total mass truncated 

system provides 

x]{t) < M!!,,Jt) < MV(0), (3.13) 

which means that all Xj, j = —n ,..., n are bounded for fixed n and not only those between 
two traps. 

Lemma 3.7 (Existence for the truncated system). For any n G N and initial data x"'(0) = 
{x'J{0))j=-n,...,n such that 0 < x”(0) < oo there is a solution to (3.12). 

Proof. Since (3.12) is an ordinary differential equation with locally Lipschitz continuous 
right hand side in (0, oo)^”’"*“^, Picard’s Theorem yields a time 0o > 0 and a unique solution 
x" G (^^([Oj^o); (0, oo)^"’"''^). By inequality (3.13), the x” are bounded, hence the solution 
can be extended either indefinitely or until a time 0 i at which one or more of the x” vanish. 
In the latter case, we have x” G (^^((O,^!)) n C'‘^([O,0i]) for all j = —n ,..., n, and the 
evolution can be continued by restriction to the living particles at time 6 i. By iteration we 
thus obtain a solution, which after at most 2n + 1 steps is identically 0. □ 

Theorem 3.8 (Local existence). Given initial data as in Assumption 3.1 and Ti < T*, 
where T* = Cd^^^ as in the persistence estimate, there exist a subsequence of (x”), not 
relabeled, and a solution x: [0,Ti) —)■ to (2.2) such that x” —>• Xj uniformly in [0, Ti) for 
any j G Z. 

Proof. Let x”, n G N be the solution from Lemma 3.7 and set x” = 0 for \j\ > n. Then 
for any fixed j G Z the sequence (x”)„,eN is uniformly bounded in T*], [0, oo)) 

and by the Arzela-Ascoli Theorem compact in C'°’“([0, T*]) for 0 < a < l/(/3 + l). Without 
relabeling we may thus extract a subsequence of n —>■ oo for j = 0, ±1, ±2,... successively 
and take a diagonal sequence to obtain x: [0, T*] —)> such that x” —)■ Xj in C'°’“([0, T*]) 
and Xj G T*]) for any j G Z. By a further subsequence and diagonal argument 

we may assume that each sequence (T”)neN converges to some 6 j G (0,oo] as n —>• oo. 

We aim to prove that x is a solution to (2.2) by passing to the limit n —)> oo in the 
integral equation 
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for Q < ti < t 2 <Ti and any j G Z. To this end, we fix 77 = rj{l3, d, L, Ti) > 0 such that the 
vanishing estimate holds with to = Ti and infer that for each j and n > \j±\ we have either 

r” > T* and x'j{t) > rj for all t G [0, Ti] (3-1-4) 

or 

rf < T* and x]{t) > C{t^ - (3.15) 

Hence, we may choose another subsequence of n —>• 00 such that for each j exactly one of 
(3.14), (3.15) holds for all large n. In the first case we trivially have x” > t/ for all large n 
and obtain 

ixj)~^X{x">o} xJ^X{xj>o} in -b^( 0 ,ri). 

In the second case we find that 6 j < T* is the unique vanishing time of xj in [0, T*], because 
0 = x^{t) —)■ Xj{t) for t > 6j and 

Xj{t) > lim C'(r”' - t)'^ = C{9j - t)?^ >0 

for t < 9j. Furthermore, we obtain convergence of {x'^)~^X{x'^>q} in T^(0, Ti) from pointwise 
convergence of X{x^>o} to X{xj>o}^ the upper bound 

x]{t)~^X{x^>o} < C{t^ - t)~^X{t<T^}, 

and the generalized Dominated Convergence Theorem. Finally, the persistence of traps and 
uniqueness of vanishing times imply that each xj has no more than 4L different neighbors 
in the time interval [0, Ti], and convergence of Ao- (xp in L^( 0 ,ri) follows. □ 

We conclude this section with our global existence result, which follows from combining 
Theorem 3.8 and the lower density estimate in Lemma 3.2. 

Theorem 3.9 (Global existence). For initial data as in Assumption 3.1 there exists a 
solution x: [ 0 ,oo) —)■ to ( 2 . 2 ), which has the following properties: 

1. for any T > 0 there is a constant C = C{f3,d, L,T) such that Xj{t) < Xj{0) + C for 
all t G [0, T] and all j G Z; 

2. for any T > 0 each Xj is Holder continuous in [0, T] with exponent l/(/3 + 1) and 
ll®illco.i/(/3+i)([o,r]) < C* = C{l3,d, L,T); 

3. if Tj < 00 then Xj{t) < C{Tj — ty/^h+i) qH i g [0,rj], where C = C'(/3); 

4 . for any T > 0 there is a constant K = K{l3, d, L, T) such that Tj < T implies Xj{t) > 

K{Tj — for all t G [0, Tj] and for any j G Z. 

Proof. Set Ti = T*/2 = Cd^~^^/2, where C = C{I3), and let x: [0,Ti) —)• £“ be a solution 
from Theorem 3.8. For ei G (0,1) to be chosen below. Lemma 3.2 provides a sufficiently 
large qi = qi{l3, L,ei) so that x(Ti) satisfies the lower density estimate in Assumption 3.1 
with R{k) replaced by Rq^k), L replaced by Li := qiL and d replaced by di := (1 — Si)d. 
Therefore, we may apply the local existence result once more to extend the solution up to 
the time T 2 = Ti + C'd^"*~^/2. Moreover, the a priori estimates are true in the whole interval 
[ 0 , r 2 ), once the constants that depend on the traps are adapted appropriately; this, of 
course, makes them time-dependent. 
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With 0<e,i<l,nGNwe may iterate the preceding argument to find 

n n ^ n—1 

dn = d\{{l-ek), Ln = L^qk, = Ti + | J] 

k=l k=l k=l 

and a solution up to any of the times T„. Choosing such that dn ^ 0 sufficiently slowly, 
we can arrange that —>• oo as n —>• oo. Finally, since any T > 0 is reached by the sequence 
(T„) after finitely many steps, the claimed estimates follow once more from adapting the 
constants of the local ones. □ 


4 Non-uniqueness 


In this section we show that solutions to (2.2) are in general not uniquely determined by 
their initial data. To this end, we consider once more truncated problems, this time with 
two slightly different sequences of initial particles that have the same limit but yield two 
different solutions. More precisely, for V G N we consider initial data G containing 
particles 

= 0 for \j\ > 3N, = 1 for — 3N < j < 0, = 1, 

and 

ag = l, = = for j = 0,1,..., iV - 1, 

where {Rfp)j=o,...,N-i, P = 1)2 are two decreasing sequences of positive numbers to be 
chosen; compare Figure 1.3. By we denote the solution to (2.2) with initial data 


+ ifj = 3iV, 

where (s^) is another sequence to be specified; in case = 0 we write . 
We will prove the following result. 

Theorem 4.1. Let 


(4.1) 


ld> Id* 


In 4 — In 3 
In 3 — In 2 


Then there are sequences {Rfp)j=Q,...,N-i, V G N, p = 
properties. 


3 

2 / 3+1 < 

“ 2 


(4.2) 


1,2 and (e'^)jveN the following 


1. For any j = 0,1,... and p = 1, 2 the sequence {Rfp) converges to some Rj^p > 0 as 
N —)• oo. 


2. There are a time T > 0, functions x: [0,T] —>■ x: [0,T] —)> and a subsequence 

of N ^ oo, not relabeled, such that for all j we have x^ —)■ Xj and x^ —)■ xj in 
C([0,r]) as N ^oo. 

3. Both, X and x solve (2.2) in [0, T) with initial data x(0) = x{0) = a, where Oj = 1 
for j < 0 and = 1 , asj+i = Rj,i, 03^+2 = Rj,2 for j > 0 , but we have x ^ x. 

Remark. Inequality (4.2) is a technical condition that we use to study in a simple way the 
phase portrait of an equation associated with two adjacent small particles. We believe that 
it can be removed or at least improved by more careful considerations; see the proof of 
Lemma 4.7 and the remark thereafter. 
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4.1 Overview of the proof 

The proof of Theorem 4.1 is rather technical, so we first outline the main ideas. To this 
end, we denote as above by and two solutions that have initial data as in (4.1) with 
= 0 and > 0 to be specified, respectively; for simplicity of notation we drop the 
upper index N where possible. Existence of x and x as well as convergence to solutions of 
(2.2) along a subsequence of —)• oo follows from the a priori estimates and compactness 

arguments in Section 3. Furthermore, since we will chose —)• 0 as —?• oo, both limit 

solutions will have the same initial data. 

Theorem 4.1 will follow from the fact that the differences 


\xj{t) - Xj{t)\ (4.3) 

are uniformly positive for some fixed index j = 3N^, some t > 0 and all sufficiently large N. 
To prove this fact, we will on the one hand derive detailed asymptotic formulas for (4.3) in 
case j = 3N^, 3{N^ + 1),..., 3N. On the other hand, we will require a precise description 
of solutions near vanishing times, because many small particles between 3Ai^* and 3N will 
have vanished at any positive time t > 0. 

For the analysis to come we require that non-adjacent small particles of x and x vanish at 
well-separated times. Moreover, the details will become simpler if we choose the sequences 
{Rf^p) such that the two adjacent small particles of one solution, say xsj+i and X3j+2 vanish 
simultaneously for each j = N — 1. In Section 4.2 we show that such a choice is 

possible, and to that aim it is necessary to study the dependence of solutions on {Rf^p)- The 
main technical ingredients here are a continuity result for vanishing times in Lemma 4.3 and 
estimates of the vanishing times in a local problem where two small particles are surrounded 
by two large ones in Lemma 4.4. The existence of {Rfp) with the desired properties then 
follows from a degree theory type argument in Proposition 4.5. 

As a consequence of Section 4.2, we obtain initial data such that for j = A^*,..., A^ — 1 
the particles x^j+p and xsj+p have finite vanishing times and Tsj+p, respectively. 

Furthermore, Tj = max{r3j+i, rs^+i, r3j+2} defines a sequence such that all small particles 
to the right of 3j have vanished at time Tj whereas all small particles to the left of 3j have 
not. The key step of our proof is to show that the most important contribution to x — x at 
time Tj-i is X3(j_i)(Tj_i) —X3(j_i)(Tj_i) and stems from X 3 j{Tj) — X 3 j{Tj). Ultimately, this 
contribution originates from the initial difference x^n — = £, is iteratively propagated 

through the solutions and amplified by the vanishing of the small particles. 

To illustrate the iterative step from Tj to Tj-i let us assume that at time Tj both x and 
X are equal except for the particle 3j and that we have the simplified situation 

x^k{Tj) = X'ikiTj) = 1, x^k+p{Tj) = x^k+p{Tj) = Rk,p ioi k < j, 

X 3 k{Tj) = X 3 k{Tj) = 1, X 3 k+p{Tj) = X 3 k+p{Tj) =0 for A: > j 

and 

X3j{Tj) = 1, X 3 j{Tj) = 1 + 6 (4.5) 

for some small 5 > 0. Then we obtain in Theorem 4.9 that the difference \x 3 j(Tj)—X 3 j{Tj)\ = 
6 is transported to the particles 3{j — 1) where we have 

|x30_i)(r,-i) - X30 _i)(T,-i)| ~ ^ (4.6) 

with Rj-i comparable to both Rj-i^p, p = 1, 2. Formula (4.6) quantifies the amplification 
effect provided that 6 is sufficiently small compared to Rj-i and we use it in Section 4.5 to 
choose initial perturbations for N > N^: such that 


6 n = £ 


N 


6j-i 


D4/3+le. 

^j-1 


r\j 


1 

3 / 3+1 


for j = Ai* -I- 1,..., A^ 
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indeed leads to an iterative amplification. More precisely, we show that we can prescribe 
> 0 sufficiently small and independent of N such that there is {£n)n>n, with the 
properties that 


e 


N 




.( 3 / 3 + 1 )^-^* 


and 

I ^3N, i'^N, ) “ i'^Nt ) I = 


for all large N. Then, taking the limit V —?■ oo yields the non-uniqueness result. 

The main part of the proof is the analysis of one iteration providing (4.6). It is carried 
out in Section 4.3, where we study the local problem of two small particles surrounded by 
a large one on each side, and in Section 4.4, where we deal with the full particle system. 
First, considering again the setting (4.4)-(4.5), we compute the asymptotic formulas 


X3{j-l}+p{t) ~ “ t) 


1 

13+1 


as t —)■ r3(j_i)_|_i 


by a straightforward local analysis in Lemma 4.6. Then, since x{Tj) and x{Tj) differ by an 
amount of order 5 it is natural to linearize 


Dp{t) = X 3 Q_i)+p(t) - X 3 Q_i)+p(t) ~ Scj)p{t) higher order terms (4.7) 

for t > Tj and p = 1, 2. This is done in Proposition 4.15, where we in particular obtain that 
(ftp becomes singular near the vanishing times of the particles 3(j — 1) -|-p and that its most 
singular contribution is of order 

3/3 

(^3(i-i)+i - t) '^+1 . (4.8) 

Hence, by (4.7) and (4.8) the linearization is valid only up to times such that 

(r3Q_i)+i - t) /3+1 ~ 6 , 

where x^j_i^^p and (j)p become comparable, and for later times we have to study the full 
nonlinear system with suitable matching conditions in order to approximate x^j_i)^p. The 
local version of the latter problem is analyzed in Lemma 4.7 and the consequences for 
the vanishing times in the particle system are given in Corollary 4.17. In the end, the 
approximations for x^(^j_i-j^p, X 3 Q_i)_|_p as well as f^j_i'^^p, r 3 (j_i)_|_p enable us to calculate 
(4.6) in Lemma 4.20, Corollary 4.21 and Proposition 4.22. 

Besides the key steps just described, Section 4.4 contains several lemmas that provide 
estimates for particles other than 3{j — 1) +p and their differences, which are not necessarily 
zero as in the simplihed setting above. The contributions of these particles, however, are 
small compared to (4.6), and their estimates usually follow from straightforward computa¬ 
tions and standard ODE arguments such as Gronwall’s inequality. 


4.2 Choice of the small particles 

With 0<7<l/3tobe hxed later we set 

Rj = +, j G N, (4.9) 

and let {R^, 2 )j=o,-,N-i satisfy 

IRj < R^p < IRj 


(4.10) 
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for j = 0,..., iV — 1, G N and p = 1,2. Moreover, we assume that (e'^)AreN is a sequence 
such that 

le'^l < - and lim = 0. (4.11) 

8 N^°° 

Then, the arguments that led to existence of solutions apply here with d = 3/4 and L = 3, 
and for later reference we summarize the corresponding results as follows. Recall that C, c 
denote generic constants that depend only on /3, if not stated otherwise. 

Proposition 4.2. There is a time T* > 0 such that for any N £ N, any sequences [Rfp) 
as in (4.10), and any (e^) as in (4.11) there exists a solution : [0, oo) —)• to (2.2) 

with initial data (4.1) and has the following properties. 

1. For all t G [0, T*] and all initially large particles, that is, k = —3N ,..., 0, k = 3j with 
j = 0,..., N — 1, and k = 3N, we have {t) > 1/2. 

2. There is a unique vanishing time G (0, r*]u{cx)} for each particle k = —3N ,..., 3N 
and there are two constants C,c such that if < oo for some \k\ < 3N we have 

(4.12) 

and 

x^(t) > c{t^ - t)^ (4.13) 

for all t G [0, r^]. 

3. For any Iq G (0,r*) there is rj = r]{/3,to) such that Xkitf) < r] for some |/c| < 3N 
implies < to + (a(fc(^o)/c)^^^. 

The existence of {Rfp) such that adjacent small particles of the solution x vanish simul¬ 
taneously follows from the continuity of vanishing times and their asymptotic behavior as 
particle sizes become locally small. We consider these issues in the following two lemmas. 

Lemma 4.3 (Continuity of vanishing times). There is N,, > 0 such that for any N > N,, 
and any solution x with (Rj), (Rfp) o.nd {e^) as in (4.9)-(4.11) we have 

< T* for j = N,,...,N-l, p = l,2. 

Moreover, the change continuously with {{Rfi^i)m=Nt,...,N-i, {Rfi, 2 )m=N,,...,N-l) G 

nf=l{hRk,iRkf. 

Proof. The first claim is an immediate consequence of (4.9), (4.10) and Proposition 4.2. To 
prove continuity, we hx and let x be a solution for initial data {af) with {Rfp) and 

as in (4.9)-(4.11). Moreover, we denote by y a solution for nearby initial data where 
Rfp is replaced by Rfp -|- 6 j^p with \6j^p\ < 5 for j = A^*,..., A^ — 1 and p = 1, 2. We write 
r and 0 for the vanishing times of x and y, respectively, and order (rsj+p) by size via 

0 < r3ji+p^ < T3 j,'2+p2 < ... , 

where (jm)m=i,...,A-A* C {A^*,..., A^ - 1} and {pm)m=i,...,N-N^ C {1,2}^“^*. 

Given small e > 0, our first step is to show by a continuation argument that the 
differences Dk{t) = yk{t) — Xk{t), k = —3N ,..., 3N remain arbitrarily small for t < r 3 j^+p^ — 
£, provided that we choose <5 sufficiently small. Indeed, Proposition 4.2(1 and 2) implies 

1 

Xk{t) > ce(5+i 
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for all |A:| < 3V and t < — e, and assuming that ykit) > /2 we compute 


yk{t) ^-Xk{t) ^ 

= /? 

d. 





<C-\Dk{t)\. 


Using the differential equation for and we find 


ll^fc(f') ||oO ^ C ^\\Dk{.^)\\OOl 


and by integration we obtain 


\Dkit)\ < <5exp - e)^ 


(4.14) 


for all \k\ < 3N and t < — e such that yk{t) > By choosing 5 sufficiently 

small, the latter inequality holds up to some positive t and the right hand side of (4.14) 
can be made arbitrarily small. Thus, for any 0 < w < we obtain 


\Dkit)\ < u) 


(4.15) 


and 

yk{t) > Xk{t) -uj> ^eTFi > c^/j+T ( 4 . 16 ) 

for all |A:| < 3V, all small 5 depending on e and w, and all t such that yk{t) > 

By (4.16) the above reasoning extends to all times t < — s. We conclude that 

(^k > x^h+pi ~ moreover we can use (4.15) and (4.12) to transfer an upper bound 

from xsji+pi to ysii+pi, namely 

1 1 

y3ii+pi('n3ii+pi - e) < a:3ii+pi('r3ji+pi - e) + u < Ce/s+i + w < Ce/^+i. 

Proposition 4.2(3) with to = T3ji+pj — e then yields 03ji+p^ < + Ce. 

Let now s G {2,..., V — W} be the index such that 


T5ji+pi — Xzj2+p2 — ■ ■ ■ — < Tojg+p^- 


Estimating Osj^+p^, m = l,...,s — las above, we can arrange that for given e > 0 we have 
\G3jm+Pm - T3jm+Pm I = \^3jm+Pm “ TOjl+p^ \<£ fol' 771 = 1, . . . , S - 1 (4.17) 

and 

roj^+Pm > T3js+Ps >T-3n+pi+V£ for m = s, ..., N - N^. 

The lower bound (4.13) and inequality (4.15) imply 

X3j^+PmiT3n+p^ - e) > and y 3 j^+p^{T 3 j,+p^ - e) > 

for m = s, N — and all sufficiently small 5. Using (4.12) for ysj^+p^ we find 

^3jm+Pm — TSji+pi + Cy/e > + e. 


and (4.13) gives 

y3j^+Pr^it) > C(03,^+p„ > C-eV(/3+i) 

for all t G [r3j4+pj — e, r3jj+p^ + e\ and m = s,..., N — N^. Consequently, we have 


/ 

Jt: 


y3j, 


^+p^{t)-^dt < = Ce^~V+r) 


"'"Sii+Pi £ 
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for m = s,..., N — Ni,, while for m = 1,..., s — 1 the bounds (4.17) and (4.12) imply 








r^3jm+P m 


C(03, 


'jm +p m 


- t) /3+1 dt < Ce^ /3+1. 


’’'Sji+Pi S 


Furthermore, the large particles k < 0 and k = 3j, j = 0,..., clearly satisfy 

f 

J T9. 


fTSii+pi+e /9 

yk{t) ^ dt <Ce < Ce^ . 

''''3ii+Pi~£ 


Using these inequalities in the equations for all particles yk that have not vanished at 
time + e, that is, for all k ^ {k = 3jm + Pm : m = 1,..., s — 1}, we obtain 

l-_S_ 

\yk{T3j^+pi + e) - yk{r3jj_+pi - e)| < Ce , 

and since a similar computation applies to each Xk, k ^ {k = 3jm + Pm :m = l,...,s — 1} 
we deduce 

l^ki^ji+pi + e)| < l-^A:('U'i+pi ~ ^)l + < io + Ce z^+i. 


Thus, at time +e we have x^j^+Pm = V^jm+Pm = 0form = l,...,s — 1, while the other 
particles are positive and arbitrarily close to each other if 5 is sufficiently small. Repeating 
the preceding arguments we obtain the asserted continuity after at most 2N steps. □ 

To study the vanishing times as particle sizes become small we consider the local problem 

ii =-2y-^ + y2"^ + Fi, ^ ^ 

« % 4.18 

Ta = -2Y-^ + Y-^ + F2, 

where Fi and F 2 are continuous functions of time. Solutions to (4.18) are understood in 
the sense that vanished particles are removed, and their existence for positive initial data 
follows as in Section 3. 


Lemma 4.4 (Vanishing times for local equation). Let {Yi,Y 2 ) be a solution to (4.18) with 
initial data 11(0) = Ai and 11(0) = A 2 , where 1/2 < Ap <2 and |Fi| + IF2I < P- Then the 
following properties hold: 


1. For all sufficiently small p depending only on (3 both vanishing times, denoted by ti 
and T 2 , are finite independently of Ap and Fp. Moreover, li + Tl < 4 for all t. 

2. Given e > 0 there exists po = po{(3,s) ~ e such that for all p < po the conditions 
Ai — A 2 = e > 0 and |Ti| + ITII < P imply T 2 < ti and Yi{t) — Y 2 {t) > e/2 for all 
t £ [0,T2]. 

3. With u = \Ai — A 2 \ and p as above we have 


max 

p=l ,2 


A: 




/3 +1 


as (n, rj) —)• (0, 0) 


independently of Fp and Ap, p = 1,2. 


Proof. Until the first vanishing time we have 


Yi+Y 2 = + Ti + T2 < -(Vi + Y 2 ) ^ + P, 
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and since Ai -h A 2 < 4 the sum Yi + I 2 decreases for all r/ < 4 ^/2. In fact, since 
rj<(Ai + A2)-^/2, we have 

yi + y2< -^(Vi + Vs)"^, 

thus the first particle vanishes at a time smaller than C{Ai + ^ 2 )^"*"^ < Thereafter, 

the same argument applied to the remaining particle shows that the second vanishing time 
is bounded by C{Ai + ^ 2 )^’'"^ < too. 

As long Yi > y 2 > 0, which by the assumption for the second claim is true initially, the 
difference Yi — 12 satisfies 

id - 12 = -3(y/ - ^ 2 "^) + Fi - F 2 > -r/, 

and we obtain 

Yi(t) - y 2 {t) >e-rit>e- ry max(ri, r 2 ). 

Hence, the second assertion follows with r]Q = e/(2max(ri, r 2 )). 

To prove the third claim, set Idnax(i) = niax(yi(t), l 2 (i)) and Ymin lt) = min(li(t), l 2 (i)) 
and denote the corresponding vanishing times by Tmax and Tmin, respectively. Letting 0 < 
u < be arbitrary and setting r = (1 — w)ldnax(0)^+^/(/3 + 1) we are going to 

show that r < Tmin < Tmax < u + 0{uj) as (u, rj) —>• 0, which yields the result due to 
|Hmax(0) - Ap\ < \ Ai - A 2 I = U for p = 1, 2. 

First, since Vnax < 4 we have < oj for all f > 0 and all rj < 4~^u}. Using 

Umax > Umin we thus find 

Umax(i) > -Umax(i)"^ “ r? > -(1 + a;)ymax(i)"^ 


for almost all t < Tmin and conclude 

Umax(t) > (Uma.(0)^+^ - (/3 + 1)(1 + u)t) ^ > ymax(0)a;?^ 

if t < min(T, Tmin)- Next, as Umin(O) = Umax(O) — 12 we have 

Umin(i) > 5 (Umax(t) “ u) > ^ (^Umax(0)a;4+T - > jUmax(0)w4+1 

for small t > 0 and all u < ldnax(o)‘^^'^*'^'''^V2- For such t we find 

( \ r 

Umi^ - yj^) + ^ = 3/3 dy + V 


< c- 


1 


Umax( 0 )/ 5 +ia ;2 
and choosing u and r/ sufficiently small we infer that 


(Unax Unin) T 1 ]: 


Umax(i) - UminW < C {u + TJU^) WoFTO ^ (4.21) 

Combining (4.19) and (4.21) we arrive at 

^^in(^) = Umax(i) + Umin(t) “ Umax(i) > Umax(0)a;/5+1 (1 -UjP+i) > |Umax(0)w4+1 , (4.22) 

which implies that by continuation the above considerations extend to min(r, Tmin)- There¬ 
fore, r < Tmin and we have 


min( T,T„ 


(4.19) 


(4.20) 


Umax (0)^+1 

13 + 1 


0{uj) + Tmin U Tjj 
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Moreover, we find 

Y~^ - Y~^ —B ^-(^+1) dv < BY~^^^^\y —Y-) 

■^min -‘max ^ ^ '^-‘min Wmax ^min; 

a Yma^-Y^ p Cuj^Y-^ 

— ^ \ \Y I V max ^ max 

\ Lmin / Lj^ax 

using (4.21)-(4.22) to estimate the first fraction and (4.20)-(4.21) for the second. From this 
and with r/l^max < a; in the equation for l^max we deduce 

Lnax(i) < -h"max(i)"^ + < “ (l - Coj^ - Y^^^{t)~>^ 

for t < T. By integration we obtain 

h"max('r)^+^ < l^max(0)^+^ - (/3 + 1)(1 - C'cai^+T _ = Y^i,^{of^^r{oj) 

where r(ti;) —)• 0 as cu —)• 0. The argument that proved the first part of the lemma but with 
r] depending on u, now yields Tmin < "Tmax <t + C'ymax(0)^+V(a;)^+h □ 

Finally, we prove that the sequences {Rfp) can be chosen such that neighboring small 
particles of the solution x for = 0 vanish at the same time. 

Proposition 4.5. There exists N^, > 0 such that for all N > N^, and any sequence {Rj) as 

in (4.9) there are two sequences {Rfp)j=o,...,N-i, P = 1,2 satisfying (4.10) with the following 
properties. 

1. The vanishing times T^jj^p, P = Ij 2 are identical for j = N^,..., N — 1. 

2. We have R^^ = Rj for j = 0,..., N — 1 and R ^2 — J = 0,..., iV* — 1 as well 

as 

^ for j = iV*,..., iV — 1, 

where (coj) is a sequence with limj^oo ojj = 0 that is independent of (Rj). 

3. There is a decreasing sequence {uij) independent of (Rj) with 0 < ujj < l/(4(/5 + 1)) 
and limj^oo <^j = 0 such that 

p/3+l 

fN _ ^ _ 

^3j+p ^ ^ ^ 

for j = N^,... ,N — I and p = 1,2. 

Proof. For N^, larger than in Lemma 4.3 to be fixed later and N > N^, we dehne a map <1> 
from the cube Q = [—1, to as follows. Given r] = G Q we solve 

(2.2) for x^ with R^^ = for j = 0,..., N — 1, R ^2 = for j = 0,..., — 1, and 

R ^2 — -^i(l + hj^j) for J = A^*,..., — 1, where {Rj) is as in (4.9) and (wj)j>o a sequence 
of positive numbers also to be fixed later. Then we set 

‘J>(^) = (T3^+l-r3^+2)f=A- 

By Lemma 4.3 the function <1> is well-dehned, finite and continuous. We aim to show that 
the topological degree deg(<l), Q, 0) is nonzero, which implies the existence of ry* G Q such 
that 4>(?y*) = 0, that is, = "^^+2 for j = A^*,..., A^ — 1. For that purpose, we use that 
the degree is invariant under the homotopy 


< 


T(A,7y) =-(1 - A)77 + A$(ry), Ag[ 0,1], ry G Q, 
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and thus 


1 = deg(-Id,Q,0) = deg(^'(0, •),Q,0) = deg(^'(l, •), Q, 0) = deg(d>,Q,0), 


provided that 0 7^ '^{ X , r ]) for all A G [0, 1] and rj G dQ . To compute 'I'(A, r/) suppose first 
that T] G dQ satisfies = — 1 for some m G {iV*,...,V — 1} and \rik\ < 1 otherwise. 
Setting 


yi (0 = 




y2{t) = 


Rm Rm 

and = 1, ^2 = 1 ~ we aim to apply Lemma 4.4 with 

-P 


Flit) = 




Rn 


F2{t) = 


^Sm+aiRm t) 




Rn 


Since ^1 — ^42 = oom and |Fi| + |F 2 | < CRm by Proposition 4.2, we need CRm < ~ 

LOm, which due to Rj = 7 -^ can be arranged for some decreasing sequence (coj) independent 
of Rj = 7 -^ and 7 < 1/3. We thus find > f^m+ 2 ^ ^{ri)m > 0 and 'I'(A, 7 )m = 

(1 — A) + 4>(?7)m > 0 for A G [0,1]. Finally, the same argument yields 'I'(A,? 7 )m < 0 if 
1]m = Tl. 

The sequences {Rfp) corresponding to rj^ G Q, whose existence is now proven, satisfy 
the second claim of the proposition by definition of <h. The last assertion is a consequence 
of Lemma 4.4(3) applied to yp{t) = /Rj for j = W,..., iV — 1 and p = 1, 2, 

which yields 


T^j+p 

1 



/3 + 1 



0 


as 


OJj 


0 , 


and finishes the proof. □ 

Remark. Exactly as in the proof above, we can apply Lemma 4.4 also to x and obtain 


^3j+p 




/3 +1 


< ujjR- 


1/3 + 1 


for j = ..., — 1 and j? = 1, 2, where and ojj are as in Proposition 4.5. 


4.3 Auxiliary results for the local problem 

In this section, we collect some auxiliary results for the local problem (4.18), which we use 
to study the propagation of the perturbation. The first two lemmas address the asymptotic 
behavior of solutions. 


Lemma 4.6 (Power law for simultaneously vanishing particles). Given Ei, T 2 G ^([O, 00 )), 
suppose that Li,Y 2 £ ^^(O,^) H (^^([Ojr]) solve (4.18) with initial data yi(0), 1 ^ 2 ( 0 ) G 
[E/2,3i?/2] for some R > 0 and vanish simultaneously at time r > 0. Then there is a 
constant C such that with 7 = ||Ei||oo + IIE 2 II 00 we have 


Yp{t) - ((/5 + l)(T-t))^+i 


< CrjQ 


t) < CR^+^ 


for all t G [0, r] and p = 1, 2. 
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Proof. The arguments for estimate (3.2) and Lemma 3.5 (with L = 3 and d = 1/2) apply 
to Yi , Y 2 and yield 

0 < c < —<C <00 (4.23) 

(r — t)/5+i 

for p = 1,2 and all t G [0, r). We define a new time s G [sq, 00 ) by 

^/3+lg-s(/3+l) =f-t, ^/3+lg-so(/3+l) ^ f 

and set 

Wr,{s) = = ie%(t), p = 1,2. 

(r-t)/3+i ^ 

Then (4.23) implies that Wp{s) G [c, C] for all s > sq and that sq is bounded from above 
and below independently of R, because r is of order Moreover, Wp satisfies 

^tTp(s) = Wp{s) + (/? + !) (-2Wp{s)-f^ + W^.p{s)-f^) + (/3 + l)Rf^e-^^Fp{t) 

for p = 1,2, and taking the difference of both equations we obtain 

^(Wi - W 2 ) = ITi - IT 2 - 3(/3 + 1) + (^0 + l)R^e-^^ (Fi - F 2 ). 

Next, multiplying by sign(VTi — W 2 ) and observing that {Wff^ — Wff^) sign(lTi — IT 2 ) < 0 
we find 

£|Wi - W 2 I > |Wi - IT 2 I - (^ + l)R^e-^f^r^ 
where we also used ||Ti||oo + IIT 2 II 00 < V- Integration then yields 

III^i - IT2 |(s 2) > (|ITi - W2\{si) - (/? + ds^ 

> (\Wi - W2\{si) - (/3 + l)i?^7?e"®i^) (4.24) 

for all S 2 > Si > sq. We conclude that 

|lTi - IT 2 |(s) < (/3 + l)R^r]e-^^ for all s > sq, (4.25) 

because otherwise the existence of si > sq such that |ITi —IT 2 |(si) > {l + i'){ld + l)R^pe~^^^ 
for some u > 0 and (4.24) imply that |lTi — IT 2 |(s) grows exponentially for s > si in 
contradiction to Wp G [c, C]. 

We now use (4.25) and Wp > c to linearize and estimate 

< C|ITi - W 2 |(s) < 

for p = 1,2, and we deduce that 

^Wp(s) = Wp{s) -i/3 + l)Wp{s)-P + R^7je-^^Fp{s) 

where Fp is bounded by a constant depending only on j3. Standard ODE methods imply 
that the only solution of VE = W — {jd + 1)IE“^ with 0<c<IT<C'is the constant 
function VTo = (/? + l)^/(/^+^), and arguing by contradiction as above we obtain 

\Wp{s)-Wo\<C\\Fp\\^pRPe-^^ 

for all s > So or, equivalently, 

^I^-(/3 + l)?^ 

(r — t)^+i 

for all t G [0, r] and p = 1,2. 


< C'||Tp||oop(r - f)/3+i 


□ 
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Lemma 4.7 (Asymptotics near vanishing). Suppose that /3 > /?* where /3* is as in (4.2). 
There are constants > 0, 5* G (—1,1), eq > 0, r/o > 0, Tq > 0 such that the following 
result is true: IfYi, Y 2 solve (4.18) (with our usual convention that 0“^ = 0 and that Yp 
remains 0 once it has vanished) in the interval (—T, 00 ), where Fp are continuous such that 
lli^lloo < r/ < r?o, r > To and and if 


(-T) - (/? + ^ 




then the vanishing times Tp of Yp satisfy 

|ri - (5* + 1)1 + |t2 - (5* - 1)1 = o{l)p^o + o(1)t^oo- 

Proof. The proof consists of four steps. In the first three steps we construct a solution for 
(4.18) and Fp = 0 which has the desired asymptotic behavior as t —)> — 00 . In the final step 
we show that we reach every solution with sufficiently small ||Tp||. 

Step 1: Backwards solution. The equation Yi = —2Yf~^ in (—1,1) with terminal data 
11(1) = 0 has the unique solution 


yi(t) = (2(/3 + i)(i-t))4+T, 


and we consider 


Y = -2Y-^ + Y-^, 


Y 2 = -2Y-^ + Y-^ 


in (— 00 , —1) with data Yi(—1) = (4(/3 + l))i/(^+i) and T 2 (—1) = 0. A local solution which 
is positive for t < —1 exists, and since Ymw = min(li, I 2 ) and Vmax = max(li, I 2 ) satisfy 


Vnin < -V 


-h 


and 


Y > —Y~d 

J- max rviQv 


we obtain 


((/3 + l)(-l-t))^+l < ymin(t) < Vmax(t) < ((/3 + 1)(3 - t))/3+1 . 


Thus, global existence follows and we find Yi >12- 

In order to study the asymptotic behavior of Yp{t) as —t becomes large, we let 


Wp{s) = = Yp{t)e-^^ 


where 


-t = e^ 


(-t)/3+l 


s > 0, 


which satisfies 


Aw — __i_ — W~^ 


(4.26) 


Subtracting the equations for ITi and IT 2 from each other and multiplying by sign(lTi — IT 2 ) 
we deduce 

T..,/ IVL 1 -IT 2 I 


i-s\Wl-W2\<-^ 


/3 + 1 


and 


|lTi(s) - 1T2 (s)| < |fTi(0) - W2{0)\e~^ = (4(/3 + 1))?+^ e~^ 
for all s G (0, 00 ). With this estimate and the lower bound 


Wp{s) > Yr 
we linearize W(~^ — Wff^ to get 


PKoj^.mmit)e /5+1 > C(-l + e")^+ie /3+1 = C(1 - e"") z^+i 


_ H^-P| < 


-h\ 


- I'Ll < Cj-^e-ra 
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and 


^sWp = - 


_^P_ ^-/3 Q 


e /3+1 


1 - e-M 

1 

for s > 0. As in the proof of Lemma 4.6 it follows that Wp = Wp — {j3 + 1 )p+^ satisfies 

\wpis)\ < C 


1 — e" 


for s > 0, and since Wp{Q) = y^(—1) — (/? + we conclude that |rcp(s)| < Ce 

for all s > 0. Thus, linearizing (4.26) around (/3 + i)i/(/^+i) we obtain 


A ( 

ds \W2j 


-1 /2/3 + 1 

/3 + 1 V -/5 


-/3 

2/3 + 1 


Wl 

W2 


+ O (rc^ + w^'j ■ 


(4.27) 


for s sufficiently large, say s > sq. The solution of (4.27) can be written as 




(so) + 



— r)0 (u;i(r)^ + W 2 {r)'^) dr. 


where 




2 


1 _ 

+ A 


+-so) 


3/3+1 

/3+1 




is the fundamental matrix for the linear part of (4.27). The contributions from +W 2 are 
of order 


and 


rs /*oo roo 

/ dr = e"" / dr - e"* / e"" dr = e-"°e-" + O (e-^*) 

J Sq J 5q J S 

r /3 + 1 ^ n ( 

e P+^ e /5+1 cIt- =-e /3+1 g /s+i + U [e ^+1 

Jo 3/3 + 1 \ ) 


with no mixed terms appearing since the eigenvectors corresponding to both powers in T 
are orthogonal. As a consequence we have 



= Ao 



e-" + Bo 



- 3/3+1 

e /3+1 +0 



+ e 



(4.28) 


for all s > So, where Ao,Bo are two constants and the contribution of order e is present 
in Wp if and only if the contribution of order e~^ is. 

Step 2: Removing the order e~^. For a given shift S'* G M let Yp{t) = Yp{t — S^) and set 


Yp{t)e /3+1 = Wp{s) = (/3 + 1) /3+1 + rrp(s) 


for t < min(0, S* — 1) and = —t. Repeating the considerations from Step 1 we find 



A 5 * 



e-" + Bs, 



- 3/3+1 

e 13+1 +0 



+ e 



(4.29) 


for some constants As ,, Bs, and s > ss, where the ss* is chosen appropriately; as above, 
the contribution of order is present if and only if the term of order e~^ is. On the 
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other hand, writing —t + 5* = + 5* = e^(l + 5*e and linearizing (1 + 5*e = 

1 + 0:5*e“* + 0(e~^^) we deduce 

Wp{s) = Yp{t)e-^^ = Yp{t - 5*)e"i^eSf = Wp(S){l + 5*e-")i^ 

= Wp(S) (l + + 0(e-2^)) 

and thus 

Wp{s) = {/3 + l)“?^5*e“'* + Wp{T) + ^ ^ ^ S^e~'‘wp{s) + 0(6“^"*). 


Formula (4.28) for Wp and the definition of s imply 



Ape-^ /1\ Bpe ^/^+i 

1 + 5*e-* Vly (1 + s^e-^)^ 



+ 0 



+ e 


-2s 


3 / 3+1 

/ 3+1 


SO that by combining the previous two equations we obtain 



((/3 + l)-^5*+^o) 



e-* + Bp 



. 3 / 3+1 
’ /3 + 1 


+ 0 



(4.30) 


as a second representation for w. With 5* = — (/3 + and B* = Bp the term of 

order e~^ in (4.30) vanishes and thus (4.29) implies that there is no contribution of order 
e“^^, either. Comparing the orders of the remaining error terms we arrive at 



= B* 



■ 3 / 3+1 

e /3+1 



for s > sst or, equivalently, 



(/3+ l)/3il(-t)/3+l 



3/3 


+ B*(-t) Z’+I 



(4.31) 


for t < — . 

Step 3: Properties of B^, and 5*. Since Yi > Y 2 , we obviously have B* > 0. Moreover, if 
B* = 0 then (4.27) implies w = 0, which contradicts the construction of Y. Thus, we have 
B* > 0. 

Next, we show 5* < 1 by comparing Y with Y{t) = [(/3 + 1)(—Suppose for 
contradiction that 5*^> 1 so that Fi(0) > 0 and 1^(0) > 0. Then we have Yi{t) > Y{t) for 
small (—t), but also Y 2 {t) > Y{t), because either 12 ( 0 ) > 0 or 

for small (—t) such that Y 2 < Yf /2. Assume now that there is a first time to < 0 such 
that Yi{tp) = Y{tp) or Y 2 {tp) = Y{tp). Then clearly l^to) ¥" Y 2 {to), because otherwise 
uniqueness for the ODE yields Yi(t) = Y 2 {t) = Y{t) for all t < 0 in contradiction to the 
asymptotics of Yp proved above. However, if ^ 2 (^ 0 ) = ^{to) < hi (to) then we find 

0 < ^(V 2 - H)(to) = Vi(to)-^ - V(to)-^ < 0, 







Mathematical analysis of a coarsening model with local interactions 33 


and in the same way the case hi (to) = ^(to) < ^2 (to) is excluded. Thus, we conclude that 
yi(t) > y(t) and h 2 (t) > y(t) for all t < 0, which again contradicts the asymptotics of Y. 

To prove 5* > —1, we consider the phase portrait of {Zi, Z 2 ) where Zp{s) = VTp((/3 + 
l)s)/(/3 + l)V(/3+i), p = 1,2. From the previous steps we know that 


lZ, = -Zr + 2Z-^-Z-f^, 
IZ2 = -Z2 + 2Z-^-Z-^, 


(4.32) 


and that the trajectory (Zi, Z 2 ){s) approaches the point (1,1) parallel to the vector (—1,1) 
as s —)> 00 . On the other hand, if 5* < 1 we have ^ 2 ( 50 ) = 0 and 0 < Zi(so) < 2^/^^^'^'^ 
where sq is the time corresponding to t = S* — 1, and a calculation shows that on the line 
which connects (1,1) to (2^/^^+^\0) the dynamics of (4.32) point to the left of that line if 
j3 > compare Figure 4.1. As this contradicts the asymptotics of {Zi,Z 2 ), we must have 
5 * > - 1 . 

Step 4- Adjusting to initial data and right hand side. For |z/| and 0 such that 

_§_ _M+i 3/3+1 , , 

\v\<{l 3 + l)f>+^T / 3+1 ^ \e P+i 


we set 

%{t-v,9) = 9-l^^Yp{9{t-v)). 

Then, Y solves the same equation as Y and with (4.31) we hnd 


^ 1 1 /1 \ 3/3/1 

y(-T; V, 9)-{13 + I)^ri5+T f J - 


u{P + l) /3+ir /3+l+0(l/^T /3+1 ^ 


+ 


3^+1 


3/3 


9 /3+1 -1 B^T / 3+1 + O u9 / 5+1 T /3+1 


3^+1 3/3 


-1 


+ O (4.33) 


3/3 

= O 


-1 


Moreover, the vanishing times of Y are ti = u + 9{S^ + 1) and r 2 = + 0(5"* — 1). 

Let now yi, y 2 be a solution of 

Yi = -2Y-^ + Fi 

in (r 2 ,Tj) with terminal data yi(Tj) = 0 and of (4.18) in (—oo,r 2 ) with y 2 (r 2 ) = 0 and 
41 (^ 2 ) given by its evolution in (T 2 ,Ti). If r]Q is sufficiently small, existence and uniqueness 
follow from the same arguments as in Step 1, and a Gronwall argument implies 

\Ypit)-Yp{t)\<C{-t)rj 

for t < + 2 , because Yp and Yp behave like a power law near their vanishing time and are 
bounded from below once they have reached a certain size depending on (3 and p. Using 
(4.33), we therefore find 

y(-r) -ip + i)p^T^ - B^T~^ =o(jjT + . 

The claim now follows from the inequality between p and T and from the fact that varying 
u and 9 sweeps a neighborhood of [{(3 + l)r]^/(^+^)(l, 1) — □ 






34 


M. Helmers, B. Niethammer and J. J. L. Velazquez 


\ \ N, \ \ rr 

"'••''''NW \ \ \ \ I j 
\ \ \ \ \ 

\ \ \ 

N \ \ \ \, \ . . 
\ \ \ { 
\ \ 1 / 

1 

’ /1 \ ’\\ 
V/ M \ \ 
,'/ f M \ \ 

/> n \ \ \ \ 

/11 \ \ \ \ \ 
y\ \ \\\\\\\ 

•\ \ \\ \\ \\ \\ 
\\\\\\\\\\ 


i I 
1 1 
I / 




1 11 / J / 

I / / / / //-^ 
/ / / ///-^^ 
/ / / / 

/ / // 7'-^^'- 

/ / /y ^ ^ 
/y' —— 


:r::::^i 

'y/t i 




NWWv.'v'v 
•.\ V N W NV 
\‘.\\\NNN 
\ \‘.\\\N\ 
\ \’.\ N N \ 
\ \ \*.\ \ N 
^ \ \ W \ N 
\ \ \ \*. \ 
\\\ \ \ \ \ ^^l 


0 


1 


2/5+1 


Figure 4.1: Phase portrait for (4.32) and /3 > /3*. Every solution (Zi,Z 2 )(s) of (4.32) approaches 
the point (1,1) as s —>■ oo either along the dashed line or the dotted line. In particular, the solution 
constructed in the proof of Lemma 4.7 runs along the dotted line. This implies S'* > — 1 because for 
S < — 1 the initial values are in the region left of the solid line connecting (1,1) and (2^/^^+^\0), 
where the dynamics of (4.32) lead to the wrong asymptotics. 


Remark. Showing that S* > — 1 is the only part of the non-uniqueness result where we use 
P > f3^. Indeed, our argument breaks down for P < (3^, because near the axis Z 2 = 0 the 
dynamics of (4.32) do not point to the left of the line connecting (1,1) and 0). 

However, formal considerations as well as the phase portrait close to (1,1) indicate that 
still only trajectories for S* > —1 reach the point (1,1) along the correct direction (—1,1), 
yet a proof seems to require a more detailed investigation of the phase portrait, which is 
beyond the scope of the paper. 

Our final auxiliary result concerns the solution of an ODE that will appear as lineariza¬ 
tion of the original equation. 

Lemma 4.8 (Solution of linear equation). For 0 < T < f let Y = (yi,y 2 ): [T,f) —>■ 
solve the equation 


d /VA /3 1/2 -A /VA 

dtV^2y /3 + lr-tV-l ‘^jyj 



(4.34) 


with initial data Y{T) = G and continuous F = {Fi,F 2 ): [T,t) —?• Then the 

following statements are true. 


1. If Fi = E 2 = 0, then 

{t) = ^t-TyY\ 

where 


$(t;r,r) 


/3 3/3 

1 {f-T\^+^ (I 1\ 1 /f-TV+' /^+1 -1\ 

A+++J (1 U + 0 ' 


2. IfY^ = ( 0 , 0 ) and Fi,F 2 are constant functions, then 

O ^ - A - T)<I>{t-,T, f)A 
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where A = (^ 1 ,^ 2 ) is the solution of 


A 2 


n 


/3- 


2 -1 
1 V-l 2 


+ 


Fi 

F 2 


3. IfY^ = (0,0), then 


\yp{t)\<c 

for p = 1,2 and all t G [T, f). 


3/3 

T — S \ / 3+1 


T — t 


{\F^{s)\ + \F2{s)\)ds 


Proof. It is easily checked that r) is a fundamental matrix for the homogeneous 

linear equation associated with (4.34) and that ^(T;T,f) is the identity matrix. Hence, 
the first claim follows immediately and the second by a direct computation. Finally, the 
variation of constants formula 


Y{t)= f ^(f; T, r)<I>(s; r, r) ^F{s)ds= [ ^{t; s,t)F{s) ds 


implies the third assertion. 


□ 


4.4 Iterative estimates 


The aim of this section is to derive the main iterative estimates for the non-uniqueness 
result. To this end, let and be solutions to the truncated problem with initial data 
with and without some G (—1/8,1/8), respectively. For simplicity of notation we drop 
the upper index N. 

We denote by Tk and fk, k = —3N ,..., 3N the vanishing times of x and x, respectively, 
and dehne 

Tj = max{r3j+i, r3j_|_2, T3j+i} 

for j = 0,. .., iV — 1; for convenience we set T^r = 0. Lemma 4.3 immediately yields 

< Tj < 

J J 

for j > Wj and by Proposition 4.5, the remark thereafter and the definition of {Rm) we 
have 

min{r3Q_i)+i, r3Q_i)+2, r3(j_i)+i} - Tj > (4.35) 

as well as 

mm{r3Q_i)+i,r3(j_i)+2,'r3Q-i)+i} ^ ^ 

rp ^ C > i 

for j > W- As a consequence, the vanishing times of 3(j — 1) + p and 3j +p are separated 
from each other and we may consider non-adjacent small particles separately. 

In the following we restrict our attention to the time interval and we refer to 

the vanishing particles 3(j — 1) + p and their neighbors 3{j — 1), 3j as the active particles; 
the others we call inactive particles. To study the change of x — x we group the inactive 
particles according to size and location and set 


IX X11 s • — 


sup 
N,<k<j-2 
p=l,2 


|^3A:+p ^3fc+p| 


R 


4+1 


\x-x\\i^j-.= sup |x 3 fc-X 3 fc|, 
M<fc<i-2 

|x-x||rj:= sup |x3fe-X3fc|. 
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Moreover, we denote by 


Dp j . p 0 , 1 , 2,3 


the differences of the active particles. 

The following theorem is the key result of this section and summarizes the propagation 
of a perturbation from the particles 3j at time Tj to the particles 3(j — 1) in the time 
interval [Tj , Tj _ i ]. 


Theorem 4.9 (Iterative estimates). Let j3 > where /3* is as in (4.2). There exist 
5* G (0,1), iV* G N, 7 * G (0,1/3), a* > 0 and C* > 0 such that for all N > V*, 6 G (0, <5*) 
and 7 G (0, 7 *) the following property is true. If for some j = V* +1,..., V — 1 the partieles 
3j satisfy 

\Ds,,{Tj)\ = 5 (4.36) 

and the other partieles are estimated by 

\DoATj)\ < 65., \D,j{T,)\ < R^+l56., \D2jiT,)\ < (4.37) 

as well as 

||x — x\\sj{Tj) < 66., ||x — x\\ij{Tj) < (5(5*, ||x — x\\rj{Tj) < 26, (4.38) 


then there is 


6 = 6 ( 6 ) G 



(^4^+1^) sAi 





such that we have 

\Ds,,.,(Tj.,)\ = \DoATj-i)\=6 

and 


IIo,j-i{Tj-i)\ < < 5 ( 5 *, 


|T> 17 _i(T,_i)| < R^+^66., 


\D2,j-i{Tj_i)\ < R^+^66. 


as well as 


^ ^11 sj —1 (Tj— 1 ) ^ 66., 


^111 (^—1) — 66 ., ||x x \\ j . j —]^( Tj —\) ^ 26 , 


provided that C.6 < 6. 


For the rest of the section we drop the index j wherever it is not necessary and write 
II • ||s. Dp and so on. Moreover, it is convenient to abbreviate 

._ -h _ J pjh ._ -0 _ ^-0 

0 ■“ ■^30-1) -^3 •“ -^37 ‘^3i • 

If not stated otherwise, we also assume N. < j < N where N. is sufficiently large so that 
all previous results apply. 

The next lemma bounds the change of particles up to time T)_i, and Lemmas 4.11-4.12 
provide a few simple estimates for the particle differences. 

Lemma 4.10 (Constancy of left inactive particles and all large particles). For all t G 
[ 0 , Tj-i] we have 

\x 3 k+p{t) - X 2 ,k+p(A\ <CA~^~’^Rk for k = N.,... ,j - 2, p = 0,1,2, 

\x 3 k{t) - X3fc(0)| < CRj-i for k = j - 1,... ,N - 1. 


In particular, choosing 7 sufficiently small, we obtain 

X3k+p(t) ^ '^Rk 

for all t G [0, Tj_i] and k = N.,..., j — 2. Similar inequalities hold for x. 
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Proof. From Lemma 4.3 and inequality (4.35) we know that for all t G [0,Tj_i] we have 

1 

X3k+p{t) > c{T3k+p - T)_i)/3+i > cRj-i 
for fc = iV*,..., j — 2 and p = 1,2. Since clearly also 

X3k{t) >l> Rj-i 

for fc = 0 ,..., j — 1, we conclude that 

\x3k+p{t) - X3fc+p(0)| < J |x3fc+p(s)| ds < CRjf.^Tj-1 < CRj-i < C'y^~^~’^Rk 

for k = N^,..., j — 2 and p = 0,1, 2. This proves the first inequality. 

For the third note that up to time Tj-i each particle X 3 k, k > j has seen at most four 
different small neighbors whose contributions to its velocity are of size 

— t)~4+T dt < Ct^ 


nT pT 

/ y~^dt<C (r - t) 
Jo Jo 


where y denotes a neighbor with vanishing time r. Since moreover large particles remain 
large throughout the evolution we deduce 

\x3k{t) — a;3fc(0)| < C (^Rj-i + Rk + Rk- 1 ^ < CRj-i 

for all t G [OjTj-i] and k > j. Finally, for k = 3(j — 1) the particles X 3 (j_i)_|_p, p = 1,2 
vanish before the time Tj-i, whereas we have x^j^j_ 2)^2 — cRj-i. Thus, combining the two 
arguments from above finishes the proof for x, and the proof for x is identical. □ 

Lemma 4.11 (Estimates for inactive particles). We have 


\x - xlUf) < C\\x - x\\s{Tj) + C^-^Wx - x\\i{Tj) + j \D^{s)\ds, 

Rj-2 Rj-2 

Ik - xUt) < C\\x - xUTj) + C{t - Tj)\\x - xUTj) + f \D^{s)\ds 

JTi 


(4.39) 


(4.40) 


and 


\x — x\\r{t) < C\\x — x\\r{Tj) + C / |Zl 3 (s)|ds 


(4.41) 


for all t G [Tj,Tj_i]. 

Proof. From the Constancy Lemma 4.10 we know that up to the time Tj-i the large and 
small particles on the left of the active ones remain close to their initial value 1 and Rk, 
respectively. In the difference of the equations for x and x we may thus linearize and find 

\^ 3 k ~ ^ 3 k\ ^ C\x3k - X3k\ < C\\x - x||z 


as well as 


-13 -0 

^3fc+pl - 


< < c\\x - x\ 


R 


■/3+1 


for A: = iV*, ..., 3{j — 2) and p = 1,2. Consequently, we obtain 
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and 


X - X o < 


C 


dt 11 -^ ^ ^+1 l\x x\\l + 


R 


C 1 


'i-2 


p/3+l 

^j-2 


^j-2 


(4.43) 


for all t G [Tj,Tj-i], where the last term in (4.43) stems from the equation for 3(j — 2) + 2. 
Applying Gronwall’s inequality to (4.42) and using t — Tj < 1 to incorporate the resulting 
exponential terms into the constants, we infer that 


\x — x\\i{t) < C\\x — x\\i{Tj) + C / ||x — x||s(s) ds. 


Similarly, we deduce from (4.43) and t — Tj < R^l^l < 


(4.44) 


|x - x||s(t) < Cllx - x||s(Tj) 4- jTi [ Ik - a;||z(s) + |i:>^(s)| ds, (4.45) 

Rj-2 


and combining (4.44), (4.45) and 


t ps 



|x — x||<j(r') dr ds = / (t — r) ||x — x||s(r) dr < / ||x — x||s(r) dr 

JTi JTi 


we arrive at 


|x - < C\\x - x||s(7(;') + i 3 +( \\^ - x\\l{Tj) 

Rj-2 


+ (77^+^/ Ik - x||s(s) ds H -l^o(®)|ds- 


Now, a Gronwall argument for t i-A |k —x|| ds proves (4.39). The inequalities (4.40) and 
(4.41) are derived analogously. □ 

Lemma 4.12 (Estimates for active large particles I). ITe have 
\Do{t)\ < C\Do{Tj)\ + C{t - Tj) (||x - xUTj) + - x||z(r,-)) 

+ cj^ K(3(i-1)) - <(30-1)) I (k d. (4.46) 


and 


mt)\ < C\D3{Tj)\ + C[t - Tj)\\x - x\\r{T,) + C l_ |x;_^( 3 .) - < 3 ^|(k ds (4.47) 


-(3j)l 


for all t G [Tj,Tj_i]. 

Proof. Similar to the proof of the previous lemma we linearize in the derivative of Dq to 
find 

®|I)o| < c\\x - xk + C\Do\ + |x;^( 3 ( - <( 3 Q-_ 1 ))|. 

Note that (T+(3(j — 1)) may take the values 3(j — 1) + p, p = 1, 2,3 in the time interval 
[Tj,Tj-i] and may be different for x and x. With Gronwall’s inequality we get 

\Do{t)\<C\Do{Tj)\+C f \\x-x\\sds + C f |<( 3 q-_i)) - <( 3 (j_i))| ds, (4.48) 

J Tj J Tj 
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and from (4.39) with \Dq \ < C\Dq\ we obtain 

J \\x - x||s(s) ds < C{t - Tj) - x\\s{Tj) + C' 7 ^+^||x - x\\i{Tj)'^ 

+ ['\Do{s)\ds. (4.49) 

Jtj 

Combining (4.49) with (4.48) and applying Gronwall’s inequality to t i—)• |Do(s)| ds 

proves (4.46). The derivation of (4.47) is similar. □ 

Next, we consider the active small particles, and our aim is to derive an approximate 
solution formula for Dp, p = 1,2. However, since the equation for Dp contains Dq or D^, 
we need precise estimates for the latter. 

Lemma 4.13 (Estimates for D^ and D^). Suppose that 

|-Dp(t)| < 5 % 0 -i)+p(t), P = l,2 (4.50) 

in some time interval \Tj,T^] where Tj < T* < min r 3 Q_i)_|_ 2 ,. Then 

we have 


\D^^{t) - D^^{T,)\ < C{t - Tj) (||x - xUT,) + C 7 ^+i||x - xll^T,-)) 


+ CRj-i\D^,{T,)\+C f (4.51) 

Jt, ^30-1)+1 - ■s 


and 


\Dl{t) - Dl{Tj)\ < C{t - Tj)\\x - x\\r{Tj) + CRj-i\D§{T^)\ +C f ( 4 . 52 ) 

JTj T3ij-l)+2 - S 


for all t G [Tj,T^]. 
Proof. In 


P-^3Q_i) ^ P [ 


--(/3+1) _ -(/3+l)\ . 

30-1) ®30-l) J ®30-l) 


we estimate the first term on the right hand side by 


di^O 


<C||x-x||, + C|Zl^|+C 


<C||x-x|U + C|Il^|+C-^ 


-h —h 

^3(j-l)+l ^3(j-l)+l 

l^ll 


■^30-l)+l 

using that (4.50) implies X 3 (j_i)_|_i > X 3 (j_i)_|_i /2 and thus 


-P _ —P 
%-l)+l ^30-l)+l 


l^ll 




30-l)+P 


4+1 


In the second term of (4.53) we use 

_ -(/3+1) _ (--P \ TC _ / /3 \ 

^30-1) ^30-1) “ ',^30-1)^ ',^30-1)^ 

where we have linearized s 1 —>■ , and 

1^30-1)1 — ^^3(j-l)+l 


4+1 


<c\dI\, 


(4.53) 
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noting that up to time Tj_i the particles and X 3 (j_ 2)+2 remain large compared to 

X 3 Q_i)_|_i due to the Constancy Lemma 4.10. As a consequence, (4.53) yields 


— 


and after rewriting this estimate as 


< Cl llx — xlL + 


l^il \Di 


* ' 


+ 




we use Gronwall’s inequality to obtain 
\D^,{t)-D^,{T,)\ 


30-l)+l ■^3{j-l)+l 


< Cl II. -*11, + 


X‘ 


,/3+l 

3(i-i)+i 


X‘ 


30-l)+l 


/3 

3(j-l)+l 


X 


< C 


\x - X L + 


\Di{s)\ 


_ mm 

X3(j_i)+i(s)/^+l ' X3Q_i)+i(s)^ 


+ 


exp 


C dr 


X3(j_i)+i(r)^ 


ds. 


Since the lower bound (4.13) implies 


ds /"* P 1 

It X3(j_i)+i(s)/^ Jt - ^^3(i-i)+i - m^+^ < CRj-u 


we conclude 


\D^,{t) - D^,{T,)\ <c£ (||x - x|U(s) + ds + Ci?,_i|Zlo^(r,)|. 

Using now (4.39) for ||x — x||s and Gronwall’s inequality for t i—)■ |C>g (s) — DQ{Tj)\ ds 

finishes the proof of (4.51). The proof of (4.52) is similar. □ 

Lemma 4.14 (Linearized equation for active small particles). Suppose that (4.50) holds up 


to time T* > Tj. Then we have 


d (Di 




dt \D2j (3 + 1 r 3 Q_i)+i -t \-l 


2 -1\ fDi 


-1 2 


T »2 

Wi 

a ;2 


+ 




+ 


+ (Ti) 

++}), 


(4.54) 


where 


I LOp I ^ C 


mmm ? m + m 


P+2 


p = l ,2 


(r3(j_i)+i - t)/3+i ('TsL-L+i - 

for all t e [Tj,T^]. 

Proof. We employ (4.50) to linearize the differences around X 3 (j_i)+p 

D^\ 


m 


A 

dt VLI 2 


-2 1 \ a^3(j_i)+A 

1 - 2/1 x~^ - x~^ / ^ 

/ a^3(j_i)+2/ 


L> 


Dn 


Dl 


which is 
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Then, replacing with the power law obtained in Lemma 4.6, that is 



and using the lower bound (4.13) to get 



proves the claim. 


□ 


Proposition 4.15 (Approximate solution for active small particles). Given (i* G (0,1) and 


6 G (0, 5 *) there are constants iV* = A^*(/3, 5*) and C = C{I3, (5*) with the following property. 
If (4.36)-(4.38) are satisfied for some j > then we have 



Moreover, (4.50) holds for all t G [T;',Tt_i] and we have 



where 



and 



with positive constants Ci, 6 * 2 , Cs, 6 * 4 . 

Proof. The Constancy Lemma 4.10 for time Tj and k = j — 1 yields X 3 (j_i)_|_p(Tj) > Rj-ij2, 
whereas (4.37) implies \Dp(Tj)\ < Rj-i/8 for p = 1, 2, all (5 < 5* < 1 and j > iV* provided 
that we choose sufficiently large. Thus, (4.50) is true for some T* G (Tj, r 3 Q_i)_|_i) and 
Lemmas 4.13 and 4.14 apply. Moreover, (4.37) also provides 



(4.55) 


for t G [rj,T*] and some constant M > 2 that will be chosen below. 

We now proceed to solve the linearized equation (4.54) separately for the different con¬ 


tributions on its right hand side and the initial data. Then, we determine the time up to 


which our approximation is valid and the leading order terms of the solution. To simplify 
the notation, we abbreviate 


r = T 3 (j_i)+i, T = Tj and R = Rj-i 
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throughout the proof. 

Step 1: Initial data Dp{T). According to Lemma 4.8(1) the contribution of the initial 
data Dp{Tj) to the solution of (4.54) with right hand side equal to zero is given by 

z'U«) = -f(*;r.r)("(gj). 

By assumption (4.37) and the formula for this is easily estimated as 


\zi^\ + < 2 


313 313 

(|I?i(T)| + \D2{T)\) < AR^+^S (5, 


T — t 


T — t 


Step 2: Right hand side {Dq{T), D^{T)) . Here we use Lemma 4.8(2) separately for 
(Z1 q(T), 0) and {0,D^{T)). Calling the solutions and respectively, we obtain 

after a straightforward computation that 


\Z^^\< CR^+^S 


f -T\ /9+1 


T — t 


and 


|Z^ 2 | < CR^+^5 


f -T\ /3+1 


T — t 


Step 3: Right hand side {DQ{t) — Dq{T), D^{t) — {T)). Using the assumptions (4.38) 

and (4.36)-(4.37) for p = 0,3 as well as the estimate (4.55), we obtain from Lemma 4.13 
that 


\D^{t) - D^{T)\ < CRSS^ + CMRf^+^6 J 


T — T\ /5+1 ds 


T — S 


T — S 


< CR66^ + CMRf^+^6 


f -T\ /3+1 


T — t 


and 


\D^{t) - (r)| < CR6 + CMR^+^6 ( 


By Lemma 4.8(3) the corresponding solution Z™ may thus be estimated as 


|Z™| < C6 




T — t 


R + 


f-T\ /3+1 


T — S 


ds 


< CR^^^d 


f -T\ /3+1 


T — t 


R + MRI^+^ 


Step 4- Right hand side Up. Plugging (4.55) into the estimate for uip from Lemma 4.14 
yields 


|wp(s)| < 


(^^j2^2(/3+i)j2 CMR^+^d /r - r\ 

+ 


(r — s) 


T — s 


(r —s)'3+i V'® 

Considering both terms separately, we deduce from Lemma 4.8(3) that the associated so¬ 
lutions Z‘^^ and Z^^ satisfy 


< CM2i?2(/3+i)j2 j {f- s)-w ds 

< CM2i?2(/3+i)j2 


f -T\ /3+1 


T — t 


1 

(r — t)4+i 
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and 


\Zf\< CMR^+^5 


< CMR^^^5 


f -T\ / 3+1 


T — t 


/;< 


T — s) + l ds 


f -T\ / 3+1 


T — t 


rP. 


Step 5: Continuation. From the previous estimates we find that {Di,D 2 ) satisfies 
l^pl < + |^p™| + \Zf\ + \Zf\ 


< R^+^5 


X -3^ 

f -T\ / 3+1 


T — t 


C(^d^ + 1 + R + MR^+^ + + |Z; 


•cjI I 
P I 


Since (5* < 1 we may choose £ N and M > 0, both depending only on (5, such that 

M 


C (^6^ + 1 +R + ^ 


for all j > N^. Moreover, we can estimate \Zp^\ similarly by i^M, where 0 < < 1/4 will 

be chosen below, provided that 


CMR <^^^6 


X -3^ 

f -T\ / 3+1 


1 


and due to r — T < CR^~^^ this inequality is true if 

'CM 


— < 


(4.56) 


T — t > 


-R^P+^6 


/ 3+1 

3 / 3+1 


(4.57) 


Hence, for such t >T also inequality (4.55) holds and together with (4.56) implies 

\Di\ + \D2\ < vif - t)Mi. 

As a consequence, also (4.50) and Steps 1-4 are valid up to the time t*_^ defined by equality 
in (4.57) if u is sufficiently small. 

Step 6: Solution formula for (Hi, D 2 ). Now all contributions to (Hi, H 2 ) except for Z^^ 
are estimated by 


CMRf^+^6 


. -3A 
f -T\ / 5+1 

f — t 


(6^ + R + MRP+^ + 1 / + . 


Therefore, by increasing iV* and decreasing n - both thus depending on <5* - we obtain 


Dp = Zp'^ + pp where \Pp\ < CMR^^^5 


The solution formula then follows from the linearization 


X -3^ 

f -T\ / 3+1 


T — t 


and Lemma 4.8(2), taking into account that the contribution of 0{6‘^) < 0((5(5*) can be 
estimated similar to Z^^ and incorporated into pp. □ 
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Corollary 4.16. In the setting of Proposition 15, we have 

\Dp{t)\ < CR^+l5 


3/3 

5+1 r f '^30-l) + l - Tj \ /3+1 


'^30-l)+l 


— t 


for p = 1,2 and 




+ C 4 


DsiTj] 


X3,(r,)5+ 


T ('^30-1)+! - Tj ) 


3/3 


'^3(j-l) + l - Tj \ /3+1 / 1 

^3(i-i)+i J 


+ O (r 3 Q_i)_,_i -t) + R-^ 6 


5+1A / '^30'-B+i 


- To- \ /3+1 


^30-l)+l 


— t 


+ \pi\ + I 52 I (4.58) 


for all t e 


Proof. Since D^{Tj) = 5 the first claim is an immediate consequence of Proposition 4.15, 
while the second in addition uses 

X3{j-i)+p{t) = (/3 + 1)^ (^30-i)+i - + O (r30-i)+i - t) , 

as proved in Lemma 4.6. □ 

Corollary 4.17. Assume that /? > /3*. In the setting of Proposition f.l5 we have either 


'^3(j 


/-i)+i “ 'r3(j-i)+i — {S* + 1)^* (1 + o(l)i{^_j_5.o), 


/3+1 


/3+1 


'l'3(j-l)+2 - 'l'3(j-l)+l — {S* - 1)^* (1 + o(l)i{^_j_5.o) 

if Ds/Tj) = +5 or the same equations with r 3 (j_i)^i and T 3 q_i )_|_2 exchanged if D‘i{Tj) = 
—5. Here 

3/3+1 (-1^ 

A ~B + T *-"4 


4/3+1 

B,(/3 + l)W 




where 5* and are f/ie constants from Lemma f.l. In particular, since |5*| < 1 we have 


X3{j-l)+2 < Ts{j-l)+l < 'TsO-L+i 
'l'3(j-l)+l < l'3(j-l)+l < T'3 (j-1)+2 

provided that V* is sufficiently large. 

Proof. In the following we abbreviate 


if 

if 


DsiTj) = +5, 
D^iTj) = -6, 


T = 


'’'3(i-l)+P’ 


T = T, 


and 


R — Rj—i, 


so that (4.58) reads 


2;3(i-i)+i 

a^30-i)+2 


1 .,1/1 


(t) - (/3 + 1)(f - *)->« (^1 j + 


4/3+1 

Ds{T) (f-r)5+r 


1 

-1 


+ O (r - t) + 


f-T\ /3+1 


T — t 


+ |pi| + \P2\ 
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Assuming that D^(T) = +5 we let 


yp{s) 


X3Q i)+p + A* 3^^+l 


and obtain 



4/3+1 

(/3 + 1)?TT(_s)]3TT Q + C4 


3/3 
(9 + 1 



_^ 3/3 

+ wi(A, (5)(-s) + a; 2 (A, (5)(-s) z^+i + /3+1 

where u!p{R,5), p = 1,2,3 tends to zero as ii —?• 0 or <5 —?■ 0. Furthermore, with the 
asymptotics 

'T-T = j^{1 + o{1)r^o), (4.59) 

which follow from Lemma 4.6, and the definition of A* we arrive at 



(/3 + l)/3+i(-s)/3+i 



3/3 

+ B*(-s)/3+l 



/3 3/3 

+ wi (i?, ( 5 ) (—s) + W2 (i?, ( 5 ) (—s) ^+3 + 0/3(12, ( 5 ) (—s) / 9 + 1 . 


Our goal is to apply Lemma 4.7 with yp, so we now check its requirements. First, (4.58) 
and the above considerations are valid for r + G [T,r*_^], which 

by (4.59) and the definition of means 


1 a _ Ji+L 

-^{R^5) 3 / 3+1 <s< 


A+ 


Moreover, we have 


where 


ds Up ~ ‘^Vp ^ + y^-p + ^p 

^3(i+p-2) • 


Ap = A+^ 


3/3 

1 r'\ /9 + 1 


By taking Al* sufficiently large we can make Fp and ujp small and arrange s G [—5, —C/A*] 
for some S that is large but of order 1. Then, the claim follows from Lemma 4.7, and the 
case (—(5) works by interchanging yi and 1/2 ■ Cl 


Corollaries 4.16 and 4.17 allow us to estimate the integrals in Lemma 4.12. 


Lemma 4.18 (Estimates for active large particles II). Let j3 > (5^,. In the setting of 
Proposition /.i5, we have 



-h _ —h 

o'+(3(i-l)) ®(t+(30-1)) 


ds < C 


and 



ds<C (^R^^+^6^ 


for all t G [rj,Tj_i]. 
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Proof. Since the proofs of both inequalities are identical, we consider only the first. For 
simplicity of notation we abbreviate 


T = 


'^3(j-l)+D 


T = r. 




and 


R = Rj—i. 


For t G [ Tj , T *] we have a^{3{j — 1)) = 3(j — 1) + 1, hence Corollary 4.16 and Lemma 
4.6 yield 




< 


cm 


X' 


,/3+l 


< 


3(i-l)+l 

By integration and using Corollary 4.17 we find 

rt 


T — t 


T — t 


3/3 




^a+{3(j-l)) ^/+{3{j-l)) 


ds < CR^+^5 (r - (r - r*)"^ 


< C 




1 

3/3+1 


For t G [r*,rj_i], however, we have to deal with the fact that (T+(3(j — 1)) changes. 
Here, the lower bound (4.13) implies 



X 


-H 

3(j-l)+P 


ds < C 



('^30-i)+p t) ds < C (r3(j_i)_,_p 



1 

/3+1 


for p = 1,2 and similar inequalities for x^j_i'j^p. Since the estimate for the possibly large 
neighbor cr+(3(j — 1)) = 3j is even better, the definition of t* and Corollary 4.17 yield 


ft 

-0 _ —0 
^a+i3U-l)) ^<7+(3(i-l)) 


ds < C 



1 

3/3+1 


for all t G [r*,Tj_i]. 


□ 


We now combine the estimates from Lemmas 4.11 and 4.12 with the consequences of 
Proposition 4.15 in order to complete the estimates for all non-vanishing particles. 

Corollary 4.19 (Estimates for non-vanishing particles). Let /? > /?*. In the setting of 
Proposition 4-15 we have 

C +C5, 

^ + C55^, 

for all t G [rj,Tj_i]. 

Proof. Using Lemma 4.18 and the assumptions of Proposition 4.15 in the estimates from 
Lemma 4.12 we find 

1 

\Doit)\ < C55, + C{t - Tj) { 55 , + + C , 

mm <CS + C{t - Tj)6 + C ^ 


|A( 0 I < 

l^3(i)| < 

X — x||s(t) < 

\x - x||i(t) < 

X — x||r(i) < 
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for all t G [Tj-i,Tj]. Then, Lemma 4.11 yields 

Ik - x||s(f) < CM* + H-^ f |Do(s)|ds 

Rj_2 ^j-2 


< CM* + C^-^55, + C* 


R 


,/3+l 

■i-2 


p/3+1 

Rj-2 




3/3+1 


k - < CM* + C{t - Tj)66^ + - Tj) 


3/3+1 


and 


Ik - ^llrk) <C6 + C{t - Tj) 

for all t G [Tj_i,Tj]. 


3/3+1 


□ 


Corollary 4.19 bounds the change of particle differences from above. However, we still 
need to study how exactly the perturbation is transported from the particles 3j to 3(j — 1) 
in the time interval [Tj,Tj_i]. To this end. Lemma 4.20 describes the transfer of mass 
from X 3 j{Tj) to X 3 Q_i)(Tj_i) and x^j{Tj) to x^j_i){Tj-i), respectively, while Corollary 
4.21 contains the corresponding result for Dq and D^. Of particular interest is the value rj 
below, which is the difference between x and x and which is responsible for the amplification 
of the perturbation. 

Lemma 4.20. Denote by 

+ H®30-l)+2('l'3(i-l)+l) */ 1'30-1)+1 < l'30-l)+2) 

“S*30-l)+l('l'3(i-l)+2) */l'30-l)+l > l'30-l)+2) (4.60) 

.0 */ l'3(i-l)+l = l'3(i-l)+2 

the scaled size of the remaining particle when the first one o/p = 1,2 vanishes. 
Then we have 

X3(j-i)iTj-i) - X3(j_i)(T,) = |x3Q_i)+i(rj) + |a:3(j_i)+2(7j') + O 

/■min(T3(j_i)+p) 


rj = 


+ 


+ 


+ 


x~^ - ^x~^ T^x~^dt 

■^3(j-2)+2 3-^3(i-l) ^ 3"^3i 


fmax(r3y_i)+p) 




3.r.-h 


/min(T 3 y_l)+p) 

rT-i 

/max(T 3 y_l)+p) 


^36-2)+2 2^30-1) ■*“ 2^3/ 


rj. h — d A- q. h +* 

■^30-2)+2 ■^3(i-l) ^ 


and 


X3j{Tj-l) - X3j{Tj) = |x3Q_i)+i(rk + |x3Q_i)+2(Tk - Vj 

/■min(r3(j_i)+p) 


+ 


+ 


+ 


3 3(i-l) 3-^3j ^-^30+1)^ 


rmax(r3y_l)+p) 

/min(T 3 y_i)+p) 

/■T-i 


1^-/5 _ 3 -/3 I -/3 Af 

2-^3{j-l) 2-^3j ^ ■^3(j+l) 


c-^ 

k(i-i) 


c-^ 

Hj 


®3L+i) 


4max(T3y_3)_,_p) 

where the maximum and minimum in the integral boundaries are taken over p = 1,2 and 
integrals with identical lower and upper boundary are 0 . 
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Proof. The assertions follow from direct computations considering the two cases < 

T 3 q_i )_|_2 and r 3 Q_i)_|_i > r 3 (j_i )^2 separately. Since the arguments in both cases are the 
same we only go through the first, where Tj < r 3 (j_i)+i < r 3 (j_i )+2 < 

For t G we have Xk{t) > 0 for A; < 3j and k = 3{j + 1), while X 3 j_|_i(t) = 

X 3 j+ 2 {t) = 0. Hence, the equations for Xk, k = 3{j — 1), 3{j — 1) + 1, 3(j — 1) + 2, 3j read 


^3(i-l) = X3(j_2)+2 ~ ^3(i-l)+l’ 

^3(i-l)+l = ^3(i-l) “ ^®3(i-l)+l ^30-l)+2’ 

®30-i)+ 2 = a^3Q_i)+i “ + x.^- , 

i3i = ^3Q-1)+2 - 2 ^ 3 / + ^3Q+1)- 

Adding twice (4.62) to (4.63) gives 

+ ^30-l)+2 = ~ 3aJ3(j_3)+3 + 3 ^ 3 /, 

and using (4.65) to eliminate in (4.61) we obtain 

^30-1) = “3^3(i-l)+l “ 3^30-1)+2 + X3l^j-2)+2 ~ 3®3(i-l) 3^3i ' 


(4.61) 

(4.62) 

(4.63) 

(4.64) 

(4.65) 


Integration from Tj to r 3 (j_i)_|_i and the vanishing of a: 3 Q_i)_|_i(T 3 (j_i)_|_i) then yield 


3^30-1)(3'3(i-l)+l) - - l^3(j-l)+liTj) + lx3y_i)+2iTj) 

, /■'^au-u+i 4-/3 1 -3 

- 3X3(j_i)+2(r30-l)+l) + / ^3(i-2)+2 " 3^3{j-l) + 3^3j (4-66) 

JTj 

Similarly, adding twice (4.63) to (4.62), using the result to eliminate (4.64), 

and integrating, we find 


®3j(33(j-l)+l) ^3j{Tj) 3a:3(j_i)_|_i(rj) + 2X3Q_i)_|_2(rj) 

- ix3Q-i)+2(r30-i)+i) + / ' 5^3(i-i) - 1 * 3 / + ^sy+i) dF (4.67) 

JTj 

Second, for r 3 Q_i)_|_i < r 3 Q_i )_|_2 and t G (r 3 (j_i)+i, T 3 q_i)_|_ 2 ) the non-vanished particles 
k = 3{j - 1), 3(j -1) + 2, 3j satisfy 

^3(i-l) = ^30-2)+2 “ ^®3(i-l) 4“ ®30-l)+2’ (4.68) 

4^3(i-l)+2 = ^3(j-l) “ ‘^^3{j-l)+2 4“ ^ 3 / ’ (4.69) 

and (4.64). Again, using (4.69) to eliminate X 3^^_^^_|_2 in (4.68) and (4.64) and integrating, 
we obtain 


3^30-1)(3'3(i-l)+2) - 3:3 (j_i)(t3Q_i)+i) = 

^ rT'3{i-i)+2 

2^30-l)+2('r3(i-l)+l) + / 

■l-Tay-ij + i 


_ 3 -/3 I + 

‘^3(i-2)+2 2■^30-1) ^ 2'^3j 


(4.70) 


and 


^3i(3'3Q-l)+2) - 3;3i(T3Q_i)+i) — 

f 

h^3{j-l)+2{'r3{j-l)+l) + 

J T- 


'f3(j-l)+2 


'r3(j-l) + l 


2 ■^30-1) 




c-^ 

'^3(i+l) 


dt. (4.71) 
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Finally, integrating the equations for and x^j over t G (r 3 (j_i)_|_ 2 , Tj-i) yields 

pTj-i 

X3ij-l)iTj-l) - X3(j_i)(r3(j_i)+2) = / ®3(i-2)+2 " + ^s/ (4-72) 

■^'^3Cj-1)+2 
rTj — i 

X3jiTj-l) - X3j(r3(j_i)+2) = / ^3{j-l) - 2 ^ 3 / + ^3U+1) d^’ 

'^'^3(j-1)+2 

and adding (4.66), (4.70), (4.72) as well as (4.67), (4.71), (4.73) proves the claim. □ 

Lemma 4.20 holds as well for x with rj = 0 identically. Taking the difference of x and 
X we obtain the following result. 

Corollary 4.21. Suppose that either T 3 (j_i)+i < f 3 (j_i)+i < T^j_i )+2 or r 3 Q_i )+2 < 
T 3 q_i)_|_i < r 3 Q_i)_|_i. Then we have 

Do{Tj-i) = Do{Tj) + ^Di{Tj) + ^D2{Tj) + rj +Ij 


and 

D3{Tj-i) = D^(Tj) + ^Di{Tj) + ^D2{Tj) — Vj + J'j, 

where rj is as in (4.60) and 

r3 
-3 

r-min(T3y_i)+p) 


T - — -I- X? -I- T? 

Xj — Xj -L Xj -L Xj 


-h 


3{j-2)+2 ^3{j-2)+2) 3^0 ^ 3^3 


-\D^^+\DUt 


+ 


rT3(j-i)+i 




T _ 3„-/3 _L 4=-/3 < 1^-0 _ Af 

•^3{j-2)+2 •^3{j-2)+2 24'3(j_i) x 3-r3(j_i) w 2 ^ 3 ^ ut 


./min(T3y_i)+p) 

+ f' ^ x~^ -x~^ - +2x~^ +-x~^-x~^dt 

^ - ■^ 30 - 2)+2 ■^ 30 - 2)+2 2 4 ' 3 (j_i) w 2 .^ 3 ^- x^j U(- 

•^'Tao-p+i 


as well as 


J,=jl+Jf + Jf 

/■min(T 3 y_l)+p) 


1 _ 4 n/3 I I 

3^0 3-^3 X I X 


-0 


^3(i+l) ®3(i+l) 


-/3 


dt 


+ 


+ 


rT'3(i-i)+i 


/min(r 3 y_l)+p) 

fTj-i 


1^-0 _ 1^-0 _ 3„-/3 , 4--/3 , -/3 _ --0 

24'3Q_i) 3‘^3(i-l) 24'3i 4^ 3-^37 ^ U.L 




1^-0 — ^~0 _ 3 -/3 I 0--/3 I -/3 _^-0 

2‘^3(i-l) ■^30-1) 2'^3i ^ 


Using our previous estimates for the various terms in Corollary 4.21 we can finally 
characterize D^iTj^i) and D^iTj^i). 

Proposition 4.22. Suppose that (3 > (5^. There is a constant a* > 0 such that in the 
setting of Proposition 4-15 we have 


rj\ = a* (1 + o( 1 )r._ 3 ^o) • 


Moreover, we have 


Do{Tj_3)-rj\ < [Rf_X^5y^^" 3^0 + 0 +6,)]+55,{l +mil) 


?/3+U 


30+T 
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and 


|Z)3(T,-i) + r,| < o(1)r,_i^o + O + <5*)) + ,5 + 

Proof. First, to determine the size of rj note that after one of the small particles X 3 (j_i)_|_p, 
p = 1, 2 has vanished, the remaining one has large neighbors only. Therefore, the power law 
from Lemma 4.6 (formally doubling the small particle) and Corollary 4.17 imply 

IC'I = g(/3 + l)'’+^ |r3(j_i)+2 - r3(j_i)+i|'3+i (^{1 + 

= l(/3 + ^ (1 + o(1)r,_,^o) . 

As A* depends only on j3 and X 3 j{Tj), and as the latter satisfies X 3 j{Tj) = 1 + 0{Rj-i) by 
the Constancy Lemma 4.10 we obtain 

\rj\ = a* (l + o(1)r^_,^o) 


for a constant a*. 

Next, we estimate the other contributions to Dq and D 3 from Corollary 4.21 separately, 
and since the proofs for both differences are similar, we consider only Dq. First, by assump¬ 
tions (4.37) we clearly have 

|L>o(T,) + lDi{Tj) + lD 2 {Tj)\ < 66, + R^^l66,. 

Furthermore, linearizing in ij we have 


|Xj|<C 


fmin(T3(j_i)^p) 


1/3+1 


R 


\x — x\ 


+ |L»o| + |^3|dt 


'j-i 


and using Corollary 4.19 we find 


|X]|<C 

< c 


fmin(T3y_i)+p) 



1 

3/3+1 


[Rf_t^6) (^1 + + C,5(l + 5,) dt 

(7"(^+i) + Rj-l) + CR^+l6{l + 6,). 


Finally, X? and X| are easily estimated using the Constancy Lemma 4.10 for X 3 (j_ 2)+2 and 
the fact that all other particles that appear in the integrands remain large. The upshot is 
that 


\Rj \ + |X|| < (CR-f^ + 1 ) (|r3(j_i)+i - T3 (j_i)+i| -h |r3(j_i)+2 - r3(j_i)+i|) 

< CRjf^ ^ 

= c(^Rf+^6)^^ [rI,6)^\ 

where we have used Corollary 4.17 in the second inequality. □ 

Theorem 4.9 is now a consequence of Corollary 4.19 and Proposition 4.22, once we 
choose 7 * sufficiently small and N,, C, sufficiently large. 
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4.5 Proof of Theorem 4.1 

We now look for an appropriate sequence {e^) and send N ^ oo. Clearly, by restricting 
ourselves to a not relabeled subsequence we may assume that R^p converges to some Rj^p > 
Rj/2 as N ^ oo for all j G N and p = 1,2. 

Suppose for the moment that we already have some for large N and that we can 
apply Theorem 4.9 iteratively to obtain a sequence 5j-i = 6 { 6 j) for j = + 1,N and 
Sn = ■ Then we have 


vd* (R 


j4/3+l^ ^ 3/3+1 


< Sj -1 < tO* ( R 




and after taking logarithms these inequalities become 

+C 1 + C 3 In Ri_i < In + C 2 + C 3 In R,_i, 


3/3 + 1 


3/3 + 1 


(4.74) 


where Ci = ln(3a*/4), C 2 = ln(5a*/4) and C 3 = (4/3 + l)/(3/3 + 1). Iterating (4.74) we 
arrive at 


Ine 


N 


{3/3 + l)^+i-3 


— + Si{N,j) < ln 6 j-i < 


Ine 


N 


(3/3 + l)^+i-3 


- + S 2 {N,j) 


(4.75) 


where 


O ^^\T O _ 'W' ^3 ^^Rj+k -1 _ Cp + C' 3 (j + /c — 1) In 7 

- 2^ (3/3 + 1)^ “ ^ “ 


(3;0 + l)^ 


k =0 k =0 

due to Rj = 'jR For fixed j the sums Sp{N,j) converge absolutely as —)• 00 and we have 

|5p(iV,j)| <C + C|ln 7 |(j + 1) (4.76) 

for all N > j with constants independent of N. In particular, for j = W + 1 we have 


Ine 


N 


j^ + Si{N,N* + l)<ln 6 N, < 


Ine 


N 


— + S 2 {N,N, + 1), 


(3/3 + l)^-M ' - (3/3 + 1)^-M 

which shows that for fixed Jtv* > 0 we can find e'^ > 0 for iV > W so that 

Ine^ ~ (3/3 + i)^“M ln<5Ar^ 


(4.77) 


and so that the iterative procedure Sj-i = S{5j) for j = W + 1 , • • •, 3V yields the prescribed 
value Jtv*) provided that the assumptions of Theorem 4.9 are satisfied. Moreover, from 
(4.75) we know that 


.N 


~ [ 2 ,_|_ i)Af+i-i ~ 'S'2(3V,i + 1) 

< In < -3/3 


lire 


N 


(3/3 + l)^+i-i 


S3{N,j + l) + S2{N,j), 


and using (4.77)-(4.76) we infer that 


5j-i 


(3/3 + l) 3 -i+M I In \ - C{j + 1 ) < In ^ < (3/3 + l) 3 -i+M 1 | + c{j + 1 ) (4.78) 

for j = W + 1 , • ■ ■, 3V — 1 . Since the first term on the left and on the right of (4.78) grows 
faster in j than the second we can choose 6 n^ so small that | ImJArJ 3> C'(W + 1 ) and 


In 


6 j-i 


> InC* 


Sj — 1 + Sj 
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for j = V* + 1,..., V — 1, where C* is from Theorem 4.9. We also obtain < 5n-i by 

using (4.77) in (4.74). Hence, for sufficiently small and chosen appropriately as in 
(4.77), Theorem 4.9 applies for all j = iV* + 1,..., V and all N > W + 1. The iteration 
over j yields 

~ > 0 (4.79) 

for all V > W + 1. 

Using our compactness results from Section 3 as well as —?■ 0 as V —)• oo, we pass 
to the limit along a subsequence of N and find two solutions x and x of (2.2) with the 
same initial data. Moreover, since cRn, < < CRn,, we also have —)> Tjv* > 0 as 

V —)■ oo, and (4.79) yields 


\x3N,iT]Sf^) - X3Ar,(rArJ| = > 0. 

Then by continuity of x^n^ and xsat, both differ in a time interval around T/v*. 
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